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ABSTRACT 


The chromatic number problem is to cetermine the minimum 
numoer of colors to assign to the vertices of a graph such 
that no connected vertices are assigned the same color. 
This oaper presents a branch and bound algorithm for the 
solution to the chromatic number problem and investigates 
five different branching functions. Additionallysr a method 
of coloring very sparse graphs is presented which divides) a 
graph into pbiconnected components and reduces the time re- 


quired to color the graph. 
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Peeler NObUC FILON 


Informally-s a graph is a collection of points, called 
nodes or verticesSs some of which may be connected by lines, 
called edges. A solution to the graph colorability problem 
specifies a color for each node such that any two nodes 
which are connected by an edge do not have the same color. 
The optimal colorability problem is to assign colorss as in 
the graph colorability problem, with the restriction. that 
the fewest nmumber of colors possible is used. 

For centuries, mathematicians have studied graphs and 
their characteristics with the awareness that many engineero 
ing and mathematical problems could be modeled by a graph. 
Then questions about the engineering problem could be 
answered by studying the characteristics of the graph and 
uSiNg knowledge derived from graph theory. 

One of the earliest and best documented graph colorabil- 
ity problems dates back to 1852, when a student of De Morgan 
named Francis Guthrier proposed that 4 colors was all that 
was needed to color a map such that no two bordering coun= 
tries are the same color(Ref. 1]. Note that this problem is 
transposed into a graph colorability problem by using a node 
for each country, and connecting 2 nodes with an edge if the 


two corresponding countries share a common border. 





In facts beginning in the 19th centurys, many af the en- 
gineering problems in areas such as schedulings communica= 
tions, transportation, electronics, and chemistry have been 
studied through the use of graph theory. More recently, the 
social, biologicals and environmental sciences have relied 
on graph theory to investigate specific problems. In the 
decade of the seventiess comouter science and graph theory 
have been closely intertwined as computer scientists realc- 
ized that many of their difficult problems can best be uno 
derstood by graph modeling. One of the most notable is the 
scheduling problem where two or more processes must compete 
for common resources. Ways to efficiently schedule these 
processes are better understood when the system is modeled 
by a graph. 

The intent of this research has been to take aie single 
problem from graph theoryr namely optimal graph colorabili- 
ty, and to investigate the use of the digital computer to 
find optimal colorings for various graphs. Though a seem 
ingly trivial problem on a computer, it is complicated by 
the fact that in the worst case, all known methods of find= 
ing solutions to problems of this type reauire exponential 
time. That iS as the size of the graph gets larger, the 
time required to find the optimal coloring grows exponen- 


tially. 





Tyoically, the researcher investigating an exponential 
time problem, whose size 1S very large, 1S faced with one of 
two choices: 1) he must accept a solution which 1s not 
necessarily optimal, but Some approximate or ‘'close' solu=- 
tions or 2) he must relax some of the restrictions on the 
problem. Most of the research done on the graph colorabili-= 
ty problem has been in the area of finding aporoximation alo 
gorithms, that 18, algorithms which attempt to find a rela- 
tively few number of colors without undue computational ef-= 
fort being expended([Refs. 2,5,4)]. The primary thrust of 
this research has been to investigate methods of finding 
only optimal colorings, and to study the effects of dif 
ferent procedures on the time and soace (computer memory) 
required by each. 

Specifically, the branch and bound technique 1s used to 
find the ootinal coloring of a large number of randomly gen= 
erated graohss and the branching function 1S varied to 
determine the relative value of different heuristic branch- 
ing functions. 

In summary, this research encompasses three related, yet 
diversified disciplines. From mathematics comes graph 
theory; from operations research comes the branch and bound 
techniques; and from comouter science comes problem complex- 


itye Background information is provided for these three too 





mes im the first part of this report. Following that is a 
discussion of the research done on the computer, including 
the algorithms used, the branching functions investigated, 
and the graph generation techniques. Finally, there is a 
discussion of the results obtained, the difficulties encoun- 
tered, and the ideas and notions derived during the course 


of the research. 
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II, BACKGROUND 


A. GRAPH THEORY 

While not attempting in any way to fully cover the field 
of graph theory, an introduction to the notions and termi= 
nology which relate to the colorability problem is necessary 
for an understanding of the discussion which follows. The 
definitions used here are standard in most of the literature 
on graph theorys so the reader with some background in the 
field can skio this section without any loss of continuity. 

Formally,s a graph G = (V,E) consists of ae finite, 
nonempty set of nodes or vertices Vy, and a finite, possibly 
empty set of vertex pairs called edges E—. For the purposes 
of this research the set of edges is restricted to not con- 
tain any elements of the form (x,x). That 318 no edge can 
exist from a vertex to itself. In the case where the edges 
have directions they are defined by an ordered pair of vere 
tices and the graph is called a DIGRAPH. This research has 
been strictly limited to non-directed graphs. 

The general graph colorability problem is formalized as 
follows: Given graph G and integer ke does there exist a 
function Rs V => I where I = {1,2,r,.eeek} Such that for atl) 
(x,y) & E, R(x) # RCy)? The optimal colorability problem is 
to find the smallest integer k such that a function R_ exe 


1StSe That smallest integer found is called the chromatic 
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number of the graph G. If R(x) = Rly) then vertices x and y 
are said to belong to the same color class. Note that the 
function R partitions the vertices of the graph into k color 
classes. 

Some additional terms relating to graphs are: 

Adjacent: Vertex x is adjacent to vertex y if an edge ex- 
ists between them, that is (x,y) @€ E. 

Adjacency matrix: A eno dimensionalsr n x nm matrix, used to 
describe any 9Qraph. The indices into the matrix are 
the vertices of the graph, and the elements are a 1 
if an edge exists between the two vertices, anda JV 
otherwise. For the araphs used in this research the 
diagonal elements are all zeroes and the matrix 1s 
reflective about the diagonal. The adjacency matrix 
iS a very commons and easy to implement method of 
reoresenting a graph tin a computer. 

Degree: The degree of a vertex is the number of edges in 
cident upon ity or the number of vertices which are 
adjacent to it. The degree of vertex x is easily 
calculated from the adjacency matrix by simply sum- 
ming the 1's in row x or column x. If the elements 
of the adjacency matrix are denoted ali,j) and the 
degree of vertex x is denoted d(x) then 


Ae = > alah) 


le 





mathe A path 1S Said to exist from vertex x to another ver= 
tex y if one can start at vertex x and follow the 
edges of the graph and reach vertex ye. Formally, 3 
path exists if there exists an ordering of some of 
the vertices of G {vj 7VorecerVy} such that v, = x » 
nie ean cm On mts) <= n=l, (verve. 4) © E-. 

Connected Graph: A graph 18 connected if at least one path 
exists from each vertex to every other vertex in V. 

Doubly connected graph: A graph is doubly connected if tor 
every vertex x and every vertex y at least e oaths 
exist from x to y and no other vertices are common 
to both oaths. This 1s also called a biconnected 
graph. 

Articulation point: A vertex v is an articulation point of 
graph Gif there exists some vertex x and some ver=- 
tex y such that every path from x to y includes vere 
tex ve. Note that a doubly connected graph has no are 
ErTEUlation DoOInts. 

Complete graph: A graph G is complete if all possible edges 
are present. That is G is complete if for all ver- 
tices x and y, (xrsy) & E. 

Clique: A clique is a subset of a graph such that” every 
vertex in the clique is adjacent to every other ver- 


tex in the clique. 
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B. BRANCH AND BOUND 

Many oroblems which deal with searching for ae solution 
or set of solutions satisfying some constraints can be 
solved using the branch and bound technique. It 1S espe- 
cially useful in solving minimization (or maximization) 
problems and has been successfully used in constrained oo- 
timization probdlems”7 since the late 1950's. Although there 
are several ways of describing branch and bound, the basic 
idea is to split the solution space (branch) and place a 
limits or lower bounds, on the optimal cost of the problem 
limited to each subset of the solution space. Those subsets 
whose optimal costs do not exceed the cost of some known 
solution are then repeatedly divided and bounded until a 
solution is found whose actual cost iS no greater than the 
lower bound on any of the subsets. This solution is thus 
the optimal one. The power of branch and bound comes’ from 
its ability to leave unexplored those sunsets which are 
KNOWN NOt to contain the optimal solution. 

In order to use the branch and bound technique the prob 
lem must have an expressable set of sotutions and a cost 
function COST() which maps the solutions into nonnegative 
Integers. Usually, the provdlem calls for finding that solu 
tion s for which COST(s) is a minimum. Sometimes all solu 


tions nay be desired whose cost is no greater than COST{s). 
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Many problems require that al! the solutions satisfy 
some set of constraints which may be divided into two 
categoriesr, explicit and implicit. If a solution is 
represented by an netuple Oxy Xo eeerX de then the explicit 
constraints are ones which determine what values” the x, 'S 
may take on. The imolicit constraints restrict the ways in 
which the x, 'S in a solution relate to each other. All solu- 
tions which satisfy the exolicit constraints are said to pe 
long to the possible solution space. The implicit cone 
Straints then determine which solutions iin the solution 
space satisfy the constraints of the problem. These solu 
tions make uo the ‘feasible’ solution space. 

For the graph colorability problem, we will use the cone 
vention that @ solution iS expressed as an netuple 
(x, 7x 


i 2 


color assigned to vertex i. The explicit constraints res= 


meee x) of integers where eect and f is the 


trict the as to integer values from i COq ah. Therefore, 


the possible solution space has a 


elements(n ways of pick>= 
ing an integer from 1 ton for each of n positions). The 
imolicit constraint in this problem is: 
(iji,3;) © E => x.# x. 

q 
The size of the feasible solution space is therefore 


also ae function of the number of edges in the graph and 


their arrangement. 
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The first step of branch and bound is to systematically 
divide the solution space into. subsets. This may be 
represented by a tree organizations, where the root of the 
tree represents all feasible solutions. The children of a 
node N then reoresent the subsets into which N can be divid- 
ed. The leafs, or terminal nodes are the particular solu 
tions which make up the set of feasible solutions. 

For the graph colorability problem, the following tree 
organization was used. The root node represents all feasi- 
ble netuples. fFhen, an arbitrary vertex 1s chosen and the 
n=tuples are divided into m subsets such that subset(k) con- 
tains all the n-tuples for which Aes Ke This 1s then done 
repeatedly, picking a new vertex at each level of the tree. 
For motation purposes, a partial solution Cx) Xn rece x de 
k < me 18 an assignment of colors to the vertices lrdvrecerke 
This k-tuple 1s then used to represent all the elements of 
a subdset. Using this notation, a portion of the tree for 
the case n = 4 18 Shown in fig. is Note that the lower 
bound on the cost of a node is the number of different 
colors in the partial solution. 

A branching strategy 18s a rule for determining how a 
solution 18 to be divided into subsets. For this problem 
the branching function is the method of choosina the next 


vertex to color at each level of the tree. 
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The tree structure for the Colorability Problem. 
Figure 1 

Once the organization is decided upon, the task of the 
problem solving method is to explore, or searcn this ‘State 
space tree’ until the optimal answer state is found. Since 
the time required by the search algorithm is a function of 
the number of nodes exnvlored in the tree, the search  stra- 
tegy 1S to explore the fewest nodes possible. One way to 
facilitate this is to ‘prune’ the tree as early inthe 
search as possible. Pruning takes place when a node is 


reached which violates the implicit constraints of the prob- 
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lem. Since no legal answer state could possibly be reached 
through that node, none of its children need be explored. 
The size of the search soace 1S thus reduced by the number 
of nodes below the pruned node. Obviously, the higher in 
the tree that oruning takes place, the greater the reduction 
in the search soace. 

The second way to reduce the search space is the essence 
of the branch and bound method. Each node has a cost asso- 
ciated with its and since the cost function 1S none= 
decreasingrs, no searching need take place below any node 
whose cost 1S greater than the cost of some known solution. 
A ‘search strategy' is a rule for choosing which unexplored 
node to explore next in the tree search. Several different 
basic search strategies exist for exploring a tree, 

A deoth first search proceeds down the tree from the 
Moot, exploring the leftmost children of each node first. 
At any point where pruning takes place, the search 1S backed 
up to the first node where it can again proceed down the 
meee in a different path. Tfopstosbottom, leftetoreright is a 
good description of this method. 

A breadth first search explores al! the children of a 
node before oroceeding down the tree. This might be called 


a leftetorrights, toop=to=bottom search. 
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A 'best-first' or ‘least cost' search strategy explores 
oniy the most promising path through the tree at any given 
time. A node is explored by assigning costs to all of its 
children and placing them on an unexplored list. fhe next 
node to be explored is the one on the unexplored list whose 
cost is currently the lowest. A oriority queue can be used 
to maintain an ordering of the unexplored nodes and at the 
point where the highest priority node is an answer state 
then the Search can be terminated because that node 
represents the vest or optimal solution to the problem. 
Nodes which violate the problem constraints may be assigned 
an arbitrarily high costs or merely not placed in the queue, 
which effectively prunes the tree. This is the search stra 
tegy employed in this research ef fort. 

This has been a general description of the branch and 
bound technique. More specific details of the tree arrange- 
ment, the search strategy, and the cost function used in 


this research will] be described in Section III. 


C. PROBLEM COMPLEXITY 

The complexity of a problem is said to be polynomial= 
time if an algorithm exists for which the number of elemen- 
tary or fundamental operations 1s limited by a polynomial of 
the length of an encoding of the problem. A problem is 


polynomial*space if an algorithm exists for which the amount 
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of elementary computer storage space needed for computation 
is never greater than some polynomial function of the encod= 
ing of the problem. An algorithm for a polynomial=time 
problem is said to be a ‘good’ algorithm and all problems 
for which a good algorithm exist belong together in a class 
called Pw. It may seem somewhat extravagant to group al} 
polynomial=time algorithms into one class since polynomials 
can be quite large. However, no matter what the coeffi- 
cients arer in the limit as n gets large, every exponential 
function dominates every polynomial function. Further, 
though unexplained, experience has shown that for many of 
the problems encountered which have a polynomial-time algo- 
rithm, the solution is bounded by a polynomial of small] de- 
gree. Some examples are ordered searching which is Olin vn), 
sorting which is O(€n*Iin€n)), and matrix multiplication, 
which is one 81), 

Another group of problems are those whose best ‘known' 
algorithm require oreater than polynomial time. Some of 
these are: 1) Algorithms requiring subexponential time, such 
as AS 2) algorithms requiring exponential time, such as 
Oe"), and 3) algorithms requiring superwexponential time, 
such as ote? ). The disadvantage of nonepolynomial=time al= 


gorithms is in the explosive growth of the maximum solution 


time as illustrated by figure 2. In this table the maximum 
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computing time 


as n goes from 10 to 60. 


13s shown for different 


comolexity 


functions 


Time Sizeon 
complexity 
function 10 20 30 4O 50 60 
n 00001 00002 .00003 £#=.00004 ~00005 00006 
second second second second second second 
a 0001 .0008 .0009 .0016 .0025 © .0036 
second second second second second second 
aS £001 008 oe £064 125 BPG 
second second second second second second 
~~ 1 eae 24.3 7 S22 13.0 
second seconds seconds minutes minutes minutes 
p 2001 1.0 17.9 f2o7 35.7 366 
second second minutes days years centuries 
gl .059 58 6.5 3855 2x108 =. 1153x1023 
second minutes years Gent. cent. cent. 
Comparison of several polynomial and exponential time 


complexity functions (Ref.5] 
Figure e 
Note that the only concern is with known algorithms = and 


worst-case situations. No conclusions are made about ex- 
pected or average performance of an algorithm. 
Problems for which there exists an algorithm which can 
Quess a solution and verify its correctness in polynomial= 
time make up the class of non=deterministic polynomial=<time 
problems, NP. For examples, graph coloring is in NP because 
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an algorithm exists which can guess a coloring and check 
whether oor not it 18 leaal tin polynomialetime. Such an ale 
gorithm needs O(n) time to arbitrarily assign a color to 
each of the n vertices. Then it needs to consider each edge 
of the graoh and check whether the e¢ vertices have the same 
color. This can be done in constant time for each edge and 
Since there may be O(n2) edges, the entire verification 
takes O0(n2). 

Note in all of this discussion the word ‘known’. There 
iS not one single problem in the class NP for which it has 
been proved that the problem is not in PlRef. 6]. There are 
no polynomialsbounded algorithms known for many oroblems in 
NP, but no lower bound which is greater than polynomial time 
has been proven for. any of these probdlems. Thus, it 1s 
still an open question as to whether the class NP properly 
contains P or whether the two classes are equal. 

There 18 obviously much Interest in trying to prove 
whether a given problem @ could or could not be solved in 
polynomialetime. For if one can prove that no polynomials 
time algorithm for Q could possibly exist then there is no 
point 1n expending the effort to find such an algorithm. 
Short of proving that a oroblem has no polynomialetime algo- 
rithm, however, there is some comfort in knowing that ones 


errricuit problem 1S Somehow related or at least as diffie 
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eult as many other problems for which no one has found a 
polynomialetine algorithm. For this discussion some more 
definitions are needed. A ‘'problem' is some general quese- 
tion to be answered which usually has one or more unspeci= 
fied parameters and some specified properties which the 
answer is required to satisfy. An ‘instance’ of a problem 
is obtained when the parameters are specified. For example, 
an instance of the colorability problem is obtained by 
specifying the vertices and edges of the graph and the 
number of colors for which a coloring is desired. 

A problem @ is said to be ‘reducible’ or ‘transform= 
able' to R if an instance of Q can be transformed in polye- 
nomial time to an instance of Re. In other words, if the 
answer to problem @ is ‘yes' if and only if the answer to R 
is ‘yest and if the input parameters to problem R can be 
determined in polynomialetime given the input parameters to 
Qr then Q 18 reducible to Re. The following schematic illus=- 


trates this definition: 


Input Polynomial Input Algorithm Output 
to Q time £0. aR for R 
transform for R and @Q 





Transformation of Problem @ to Problem R 
Figure 3 
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Clearlys if Q is reducible to R and R is in the class Pr, 
then @ is also in Pe. Note that nothing is said about the 
existence of an algorithm for Rs only that ‘if!’ an algorithm 
for R exists which runs in oolynomial=timers, then one also 
exists for Q@. 

If every problem in the class NP is reducible to some 
oroblem Q then @ is said to be NP*hard. In other words, & 
is at least as hard as every problem in NP, Further, if Q 
is also in the class NP then @ is NPecomplete. The sigqnifi- 
cance of a oroblem being NPecomplete is that i f 3 
polynomial-time algorithm is ever found for any NP-complete 
oroblem then it will be known that a polynomialetime algo- 
rithm exists for all NP prodlems and that indeed NP = P. 
Further, if any problem in NP is ever shown to require 
super=p0lynomial time complexity then every NPe=complete 
oroolem must also require greater than polynomial time, 

To show that a particular problem Q is NP#=complete, 
then, one must show that it is indeed in NP and then either 
show that every NP problem could be reduced to QQ or that 
some already proven NPecomplete problem is reducible to @Q. 
In 1971, Steven Cook(Ref. 7] laid the ground work for the 
current theory of NPecompleteness in a paper entitled “The 
Complexity of theorem proving Procedures.” One of his most 


Significant results was the proof that the 'Satisfiability' 
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oroblem from boolean logic is NP=complete, and it holds” the 
distinction of being the first NP=complete problem. Referee 
ence 5 contains an excellent description of the Satisfiabil= 
ity problem as well as a proof of Cook's theorem. Cook also 
proved that a variation of the Satisfiability problem, 
called 3-Satisfiability,s was NP-complete by transforming Sa- 
tisfiability into S-Satisfiability. Then, in 1972, 
Karp(Ref. 8] oroved that many decision problems were also 
NP-complete. One of these proofs was that the 3- 
Satisfiability oroblem was reducible to the colorability 
problem. Thus the colorability problem is also NP#complete. 

Inasmuch as most theorists are in agreement that the 
NP-complete oroblems are probably super=polynomials much 
work has been done in recent years to find relaxations to 
the problems which allow them to be solved in polynomial 
time. As the graph colorability problem has many apolica@ 
tions in real@life problems, many methods of finding aoprox> 
imate solutions have been tried. In many cases where a 
graph can be used to model a real=life problem, the absolute 
minimum number of colors may not be needed, but rather, some 
close approximation might be satisfactory. Unfortunately, 
even this has oroved to be a difficult problem. Many algo- 
rithms have been devised(Refs. 2,3,4] to find approximate 


colorings, but for every one, it is possible to construct a 
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graph such that the approximation algorithm will find a 
coloring which uses many nore colors than the chromatic 
number. Though the algorithm may indeed run in polynomial 
time, the closeness of the approximation to the actual 
chromatic number of the graph may be suspect. In fact, in 
1976, Garey and Johnsonl(Ref. 14] proved the following: 

Given graoh G_ with chromatic number X(G), and 
polynomialetime approximation algorithm A», which computes a 
coloring of G using A(G) colors, there can be no A such that 
for all G: 

ACG) 
— <r fon G ee <e . 
X(G) 

In other wordss any approximation algorithm wil] always 
have some inout graph for which it can find no coloring 
which uses less than twice as many colors as actually need- 
ed. 

The primary reference for the material in this section 
1S Garey and Johnson(Ref. 5). It is highly recommended for 
the reader who 18 interested in other NP=complete problems 


and further discussion on this topic. 
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T1IT. THE ALGORITHMS 


A. GRAPH GENERATION 

One of the earliest problems faced in this research was 
that of deciding what characteristics were desirable in the 
graphs to be colored and how to generate them. As the goa! 
of the research was to contrast the efficiencies of various 
heuristics for optimal colorings, the graphs used needed to 
ber, in some sense, typicals or representative of graphs 
which occur in realelife problems. Though the question of 
what iS a tyoical graph was left largely unanswered, three 
possibilities for graph generation were considered. 

One method considered was to use hand generated graphs 
with certain built-in specific features. For example, 
graphs could be designed with particular distributions) of 
degree of the vertices. Another method would be to generate 
graphs with a given chromatic number but with varying 
numbers of vertices or edges. The major drawbacks to this 
method were the necessarily small sample of graphs and the 
failure of the graphs to be typical in any sense. This 
method was therefore not used for any of the experiments 
done but was used extensively for program testing and debug= 
QING. 

Another possibility, quickly discarded, was to survey 


some of the disciplines where graph colorability models the 
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oroblem and to use some actual graphs which occur 11n Pfreale 
life. Though this might generate some interesting results 
for those specific applications,s, it would also lead to 
results which were heavily flavored by the problem being 
considered. What was needed was a selection of graphs which 
were, im some sense, representative of the set of al] 
graphs. 

The most ae method of graph generation was to use 
the computer and to create a large population of random 
graphs on which to test the various heuristics. There is 
some benefit to using random graphs since one generally 
feels comfortable that these random graphs are fairly 
reoresentative of the set of all graphs. 

The orimary drawback to random generation of graphs is 
that the generator wil! tend to generate more of the types 
of graph which have many isomorphisms than the type with few 
isomorphisms. A graph G is jiSomorpnic to graoh H if the 
vertex numbers of G Can be permuted in such a way as to make 
the resulting adjacency matrix of G identical to the adja- 
cency matrix of H. Fig 4 shows 2 of the isomorphic graphs 
for each of the only e distinct graphs possible for 4 vers= 
tices and 4 edges. 

Since graph coloring algorithms can be influenced by the 


different types of graphs, a random generator may bias the 
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investigation of the relative performances of different 
coloring programs. fhe generation technique could be amend= 
ed by somehow reducing the probability of using a particular 
graph by dividing the probability by the number of the 
graph's isomorphisms. But determining the numoer of iso- 
morphic configurations of a given graph is not a Straight= 


forward problem and beyond the scope of this research. 


UM Uk 


Examoles of the isomorohs of 2 distinct graphs. 
Figure 4 


Another technique might be to record each graph and to 
test each mew one generated against all previous to deter- 
mine whether it iS an isomorphism of an already used graph. 
But the problem of just testing two graphs against each oth 
er 1s probably NPecomplete(though considered an open queso 
tion at this time) (Ref. 5). 

In order to get on with the research goals at hand, Vt 
was decided to accept a random graph generator which gen- 
erates graph type G more than tyoe H if G has more isomor- 


Dphisms than HH, This 318 an admitted shortfal) in this 
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research effort and it 18 hoped that the algorithms experic 
mented with are not discriminated against by the shortcoming 
of this random graph generator. 

Inasmuch as unconnected graphs are no more difficult to 
color than their largest connected component, it was decided 
to run the initial tests on only connected graphs. Two 
methods of generating connected graphs were considered: 1) 
To randomly add edges until the graph was connected and 2) 
to create a spanning tree of n vertices and n=-1 edges(the 
minimum to connect the n vertices) and then add additional 
edges at random. Since method one would require testing for 
connectivity after each edge and also would fail to generate 
very many sparse graphs(found often in real oroblems) method 


2 was chosen and implimented using the following algorithm: 


—_— 


Algorithm GENGRAPH 
1. Let VC be a set of connected vertices, initially empty 
Let VU be a set of unconnected vertices, initially 
VU = {lpepeeern 
Ce. Choose a random vertex a 
3. Move a from VU to VC 
4. While VU iS not empty 
choose a random vertex a from VC 
choose a random vertex b from VU 
record an edge from ato b 
move b from VU to VC 
5S. While Number edges < Number desired 
choose random vertex a 
choose random vertex }b 
While edge (arb) exists or if a=ob 
increment a until a = n then increment b and set a = 
if a = mn and b =n then set a= 1 andb = i 
record edge (a,b) 
End 
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At a later point in the research, it was decided that 
for truly valid test results, the graoh generator must gen- 
erate purely random graphs with no requirement for connec- 
tivitye It was discovered that the connectivity requirement 
changed the characteristics of the sparse graphs and non- 
reoresentative graphs were being used. Thereforer a much 
simplified generator was used which ignored steps 1 through 
4 of GENGRAPH and merely added random edges until the re- 
guired number was oresent. Another algorithm called DIVIDE 
(to be described in a later subsection) was then called to 
divide the graph into its connected components. 

The random number generator used is identical to Grogono 
(Ref. 9) with a modification to return integers in any range 
desired. The listing for both GENGRAPH and the Random 


Number generator 1s enclosed in appendix A, 
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8. GRAPH OIVIDING 

It has already been shown that an unconnected graph may 
be colored by dividing it into its connected components and 
Separately coloring each component. A further simplifica/ 
tion can be made, however, if even the connected components 
are divided at their articulation points, if any exist. Ree 
call that an articulation point iS a vertex which is in 
every path from some vertex x to some other vertex y. If the 
articulation point is divided into 2 vertices such that the 
path from x to y 18 cuts then the graph 1s divided into 
disconnected subgraphs. For illustration purposes consider 


mags. > and 6. 





A Graph with articulation points. 
Figure 5 
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The articulation points are vertices d and f. The 


biconnected components are shown in fig. 6. 


A Graph Divided at its articulation points. 
Figure o 
In communications and transportation problems the iden- 
tification of articulation points is very important because 
these are the ‘choke’ points, or places where the network 
becomes most vulnerable to failures. Since dividing a grapnn 
at its articulation points results in biconnected components 
which can be colored separately,s, it is also important to 
discover the articulation points in the colorabdility prob- 


lem. The number of colors needed will be exactly tne same 
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as the number needed to color the component which uses” the 
most colors. If the actual color assignments are needed, 
then the graph can be ‘pieced’ back together one component 
at a time. As each component 18 added, the color classes 
are permuted so that the articulation vertex 18S 1n the same 
class in both components of the graph. 

An algorithm to find the articulation points and the 
biconnected components of a graph utilizes a depth first 
search of the vertices 1n the graph and makes use of the 
fact that v 18 an articulation point 1f and only if there 
are two vertices x and y such that every path from x to_ y 
includes v. 

The idea of the depth first search is to visit al) the 
vertices of the graoh by starting at some arbitrary start 
vertex Ss and recursively visiting every other vertex. Given 
that the search 1S at some vertex a, it follows an edge 
(arb) to vertex be. If b has already been visited, the 
search returns to aandtries another edge. If b has not 
been visiteds, then the same method is applied recursively to 
be At the point where al! the edges incident on a have been 
examined, the search backs up along the edge that took it to 
ae The search terminates at the point where it tries to 
back uo from the start vertex. If the graph was not con- 


nected, then some new start vertex is chosen from the un 
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visited ones and the procedure repeated until all vertices 
have been visited. The edges which take the Search from one 
vertex to a vertex which has previously been visited are 
called ‘back edges.' Edges which go to an unvisited vertex 
are called ‘forward edges.’ These definitions will be need= 
ed for the algorithm description. 

Figure 7 shows a graoh with 5 vertices and 7 edges. Lf 
the depth first search is started at vertex a then the solid 
arrows represent a possible set of forward edges and the 
dashed arrows are the back edges. The labels on each edge 


indicate the order in which they are traversed. 





—_m_ = =P oer OS Se le 


A Depth first search of a Graph. 
Figure 7 
The basic idea of using a depth first search to find the 
articulation points of a graph can be seen by studying fig. 
8. This schematic represents a graph with 3 piconnected 
components labeled A, Br, and C which are joined at articula= 


tion vertices v and w. S 18 an arbitrary start vertex in A. 
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A Schematic of a graoh with 3 biconnected components and 
2 articulation points. 
Figure 8 

If the depth first search is started at vertex S in A, 
it will eventually enter into C through vertex v. By the 
deoth first nature of the searches all the edges in C will be 
traversed before the search backs up through v. If the 
edges are placed on a stack as they are traversed, then when 
the search gets back to ve all the edges in C will be on top 
of the stack. 

In order to recognize the articulation points, note that 
no vertex in C will have a back edge to any vertex which was 
visited prior to ve Sor if the vertices are numbered as 


they are visited and eackK vertex is tagged with an index 
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equal to the smallest number of any vertex which can be 
reached through any number of forward edges and one back 
edge, then the articulation points can be found. When the 
search backs up from b to a along edge (a,b), if the index 
tagged to b is not less than the number assigned ar, then a 
is an articulation pointflor the root of the tree). In other 
wordsr, no vertex below a in the tree has a back edge to a 
vertex visited before a. The following algorithm by 
Baase({Ref. 6) was used to find both the connected components 


and the biconnected components of a graph: 


Algorithm DIVIDE 


tack of vertices 

tack of edges 

top = too element of SV 

Number = array to number the vertices in the order visited 
Back = array to record the last vertex reachable 

through a back edge. 


ww 
<= 
a6 tt 
on Ww 


Algorithms 
le initialize Number(i) = 0; 
2. choose arbitrary vertex §S 
3. Number(S) = } 
4. Num = 2 // next number 
Bemestack S on SV (top = 8S) 
6. while there exists an unprocessed edge from top to v do 
stack edge (too,v) on SE 
if number(v) > 0 then backltoo) = minl(back(top) ,back(v)) 
elise do 
number(b) = num 
numt+? 
back(v) = number(v) 
stack v on Sv 
end 
end 
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7. if there is more than one vertex on SV 
let v = second from top 
if back(top) >= Number(v) then 
let edge (arb) be the top edge on SE 
while Number(a) >= Number(v) and Number(o) >= 
Number (v) 
pop SE 
else Back(v) = min(Back(v),Back(topn)) 
pop SV 


goto 6 END 


The code used to impliment this algorithm is enclosed in 


appendix A. 


C. GRAPH COLORING ALGORITHM 

The essential idea of the basic coloring algorithm is to 
build the search tree until the optimal coloring is found. 
Each node in the tree represents the placement of a vertex 
into a color class and the path from the root to the node 
reoresents a oartial coloring of the graph. The cost of a 
node is the number of different colors used for that partial] 
coloring. As the nodes of the tree are generated, by itera~ 
tively expanding from the start node, they are placed in a 
priority queue such that the first node in the queue is the 


node which represents the least numoer of colors used and is 
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deeper in the tree than any other node with the Same number 
of colors. A node is expanded by generating all of its pos- 
sible children and the node to be expanded next is always 
the first node in the queue. The stopping point for the al-= 
gorithm is thus the point at which the highest priority node 
mea 6€6ftUl!)6Ucoloring) 60of the graph. At that coint no other 
node in the tree could possibly be expanded in any way which 
results in a coloring using fewer colors. The chromatic 
number of the graph is then known to be the number of colors 
used in that coloring. Note that if all possible colorings 
are desired which use the optimal number of colors, then the 
algorithm need only save the chromatic number and to contin= 
ue expanding nodes until the cost of the highest priority 
node is greater than the chromatic number. At each point 
where a node is generated which is a full coloring, it can 
be printed out and all optimal coloring assignments will 
thus be generated. 

There are 2 primary data structures used in this algo- 
rithm. These are the TREE and QUEUE data structures. The 
elements of TREE are trees made up of nodes and the _ links 
between them which organize the nodes into a tree structure. 
A node N represents a color assignment to one vertex of the 
graph. The path from N to the root of the tree represents a 


partial coloring of the graph and the cost assigned to node 
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mers the number of colors used in that partial coloring. 
Each node has a pointer to itsS parent In the tree IN order 
to maintain the links between nodes. There are ec operations 
defined for TREE. The operation CREATE generates ae tree 
with only one node called the root. The operation EXPAND 
takes as input a tree T and a particular leaf node N of fF 
and returns a tree — which is the original tree T plus all 
the oossible children of node N. 

The nodes of the tree are imolimented using an array of 
records where each record contains the following 5 fields: 
1) the depth, or level in the tree at which it 1s located, 
2) the cost assiqned to the node, 3) the vertex number which 
was colored to create the node, 4) the color assigned to 
that vertex, and 5) a pointer to tne node's parent in the 
tree. Fhe operations on the TREE are implimented in the ale 
gorithm COLORALL and EXPAND as follows: 

COLORALL (Color all vertices of the graoh) 
1. Get a node N 

ee Get number of first vertex =v. 

3. Let N be the root node 

Nevertex Vv 

N.ecost 1 

Nedepth 1 

N.ecolor ] 

Neparent pointer = NULL 
G4. While N.edepth < n do 

expand N (see algorithm EXPAND below) 


N = first node in priority queue 
S.- Print out statistics,coloringsand chromatic number. 
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The algorithm to expand a node is as follows: 


EXPAND(N) (N 18 the node to be expanded) 

1. Choose a vertex v to be added to the partial coloring. 
(described in subsection D.) 

ee et C = {1l,2seeeeN.cost}? be the set of colors used 

in the path to node N. 

3. Trace the parent pointers from N to the root. For 
each node M in the path to the root, if an edge exists 
between vertex v and M.vertex then remove color M.color 

from C. 

G. For each color c. left in C create a new node N. and 

out it In the priority queue. 


N. evertex = Vv 
N, cost = N.cost 
N. depth = Nedepth + 1 
N. scolor = 6x 
N, -parote = N 
5. Add one new color and create one more node N.. 
N.evertex = v J 
hme .cest = N.cost +1 
NJ depth = Nedepth +! 
NY .color = Necost + ! 
eesrote = N 
END 


In terms of graph coloring, step 4 creates a child 
node for each of the ways in which the new vertex can be ad- 
ded to one of the existing color classes. Steo 5 generates 
@ node which represents putting the new vertex in a new 
color class by itself. 

The QUEUE data structure is used to organize the leaf 
nodes of a tree according to their priority. There are 2 
components which determine the oriority of a node. The 
first 1S, of course, its cost as defined earlier. The 


second component 18S used to determine the highest priority 
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between 2 nodes with the same cost. When more than one node 
is available which have the same cost then the node which 1s 
deeoest in the tree has a higher priority since it is closer 
to a complete coloring. [There are two operations defined 
for the QUEUE called ADDNODE and REMOVETOP. ADDNODE takes a 
queue Q and a node N and returns a queue Q with N in its 
proper priority location. REMOVETOP takes a queue @ and re- 
turns the node N with highest priority. 

The priority queue 1S implimented as a heap, as 
described in ref. 10. The nodes are organized into a full 
binary tree such that the root is the node with highest 
priority and every node in the tree has a higher priority 
than either of its children. When adding 8&8 node to the 
tree, it 18 put in the first available location and then 
‘sifted up’ until its parent has a higher priority. When 
the root node iS removed, the last node is moved to the root 
positions then ‘sifted down’ until both its children have a 
lower priority. The binary tree 1S maintained by an array 
of pointers q() where q(1) points to the node at the root of 
the binary tree and the children of qi) are located at 
qa(2i) and afei + 1). 

In any tree search algorithm it is important to ensure 
that isomorphic paths are not generated. For the problem of 


graph colorabilitysr two colorings are isomorphic if the vere 
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tices of the graph are partitioned into color classes in the 
Same way in both colorings. Therefore, if one coloring uses 
k color classes, then there are k! possible isomorohs of 
that coloring. In order to prevent the tree from containing 
all of these unneeded paths, one must either check for pos= 
sible isomorphs with the addition of each new node or 
orevent isomorphic path generation from occuring by the 
design of the algorithm. 

There are two essential ingredients of algorithms 
COLORALL and EXPAND which ensure no isomorphs: 1) When a 
node N is expanded, at most 1 new color is added to the 
available color set, and 2) In expanding a node Ne only i 
vertex of the uncolored vertices can be added to the partial 
coloring and that vertex must be the same for all the chile 


dren of N. 


THEOREM: Any tree search graph coloring algorithm with the 


above e oroperties will not generate isomorphic colorings. 


Proof: Assume the algorithm generates a tree with 2 iso= 
morphic paths P and Q. Then, by the definition of isomorph- 
ic colorings, for each color class in P there exists a color 
class in @ with exactly the same elements. Let vertex v be 
the first vertex colored in P and @ such that the color of v 


in P is different from the color of v in @. Call the class 


G3 





containing v in path P, on and the class containing v in GQ, 


C3 ‘ If the suboscript indicates the colors then i #4 j. 
Since only one new color could be added when v was colored, 


then at least one of a or C3 contains a vertex w which was 
colored prior to ve. Without loss of generality, let that 
class be Ge. In other words, v and w are in the same color 
class in path P. But, since v was colored differently in 
path Q and wis the same in both P and GQ, there can be no 
color class in @ containing both v and we. Therefore, P and 


Q are not isomorphic colorings. Figure 9 illustrates this 


proof. 
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Two Paths P and Q which are not isomorphic 
Figure 9 
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DOD. THE VERTEX ORDERING FUNCTIONS 

The purpose of this research was to investigate dif- 
ferent branching functions In a branch and dound approach to 
solving the chromatic number problem. For the particular 
branch and bound algorithm useds the branching function is 
the order in which the vertices are selected to be added to 
the partial coloring at each new level of the tree. In this 
research, five different methods of ordering the vertices 
were investigated. Four of the methods were static orders 
ings in which the vertices were preordered and then the 
graph coloring algorithm was called to search the entire 
tree. In every path through the tree in these methods” the 
order of choosing the next vertex was constant. In the 
fifth method, the vertex chosen to be colored next was 
dependent on the particular color assignments made prior to 
the node to be expanded. This method of dynamic ordering 
results in alarger time expenditure at each node since the 
order must be recomputed at each expansion of the tree. 

Note that the purpose of any of these orderings is to 
make the search tree as small as possible, thus requiring 
fewer nodes and fewer expansions. If the amount of work 
done to expand a node (create its children) is constant for 
each methods, then the resulting size of the search tree is a 


direct measure of the benefit of any procedure. Since the 
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height of the tree will always be the number of vertices in 
the graph, the only way to-decrease the number of nodes 1s 
to decrease the branching factor (number of children of each 
node). The optimal branching factor i8r of course, 1. With 
the algorithm COLORALL, this 1S reached when each vertex to 
be colored requires a new color. It iS approached most 
closely when the vertices are ordered such that the next 
vertex is always the one with the fewest colors available a 
be colored. Three of the static orderings attempt to accom- 
olish this ordering through the use of some heuristic func= 
tion to order the vertices. The dynamic ordering actually 
determines, at each nodes which of the uncolored vertices 
has the fewest colors available. 

Ooviouslys any graph which is complete(a clique), will 
be colored with a branching factor of ones no matter what 
order is used. The same applies to a graph which is just a 
forest of vertices(no edges). The difficulty lies in the 
midrange. Another way to look at it is that coloring is 
easiest when there are few restrictions(edges) or very many 
restrictions to the colors which can be used. Each of the 
ordering methods, therefore, really attempts to find the 
vertices with the most restrictions and color them first. 


Thus, the tree has the smallest branching factor at the top 
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and more of the tree is oruned off resulting in fewer nodes 
to search. 

Following is a description of each ordering. The list= 
ing for all orderings is in the program ORDERNODE in appen=- 
dix A. The case numbers in the switch statement correspond 
to the numbers shown here. 

1. Random 

This method of selecting vertices was used as 4a 
baseline against which to compare any benefit of the other 
methods. Inasmuch as the graph was generated randomly, no 
real work needed to be done. The vertices were merely 
chosen in the order of increasing vertex number. 

2. Degree 

This method, used in the Welsh=Powell approximation 
algorithm({Ref. ll), orders the vertices by decreasing values 
of the degree of the vertices. Recali that the degree of a 
vertex 1$ the number of edges incident upon it and is easily 
calculated by summing the elements of the rows of the adja-~ 
cency matrix. The idea behind this method is that the ver- 
tex with the highest degree has more restrictions on it when 
choosing a color so it should be colored first. That color 
used is then removed from the colors available for every 


vertex connected to it. 
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3. Dvec 3 

DVEC 3 is the name given by this author to a method 
devised by M.R.WilliamslRef. 4]. In the prior method of 
ordering the vertices by degree, the underlying premise is 
that if a vertex 1 is connected to p other vertices, then it 
will be more difficult to color than vertex j which is only 
connected to qvertices (p < gq). What that method fails to 
consider is the degree, or difficulty, of the vertices to 
which i and j are connected. For example, if vertex j is in 
a clique of size 4 and i is the center of a star with 5 
points, aS Shown in fige 1046 then certainly j should be 
colored before ir as well as all the other vertices in the 


clique. 





Example of Shortcoming of Degree method 
Figure 10 


Williams recognized that the vector of degrees of the vere 
tices was the first step towards the development of the dom- 


inant eigenvector of the adjacency matrix. Recall that 
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In the standard technique, dq; iS the eigenvector 
corresponding to the largest eigenvalue of matrix A and is 


found by iteration using: 


n al n 
dj = 2. a45* 9; 
J 


ef da is the vector of degrees, then this scheme 
does take into account the degrees of the vertices to which 
a vertex is connected. And, in fact, the number of itera- 
tions used determines how far away from vertex 1 to look in 
calculating the relative difficulty of coloring vertex 1. 

Williams also determined that the relative order of 
after a few iterations. This fact was borne out by results 
of this research. Although improvement was seen out to e or 
3 iterations, little or no changes in the size of the search 
tree was found after 3 iterations. The only effect was the 
increased computation to calculate the vectors. For this 
reason, 3 iterations was chosen as a good number to use to 
compare this method against the others. 

4. Adjacency 

The rationale behind this ordering method is that 

the vertex which is adjacent to the most previously colored 


vertices wi!] probably have the fewest colors available. 
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This 3S not necessarily truer however, as a vertex may be 
adjacent to many vertices which are all in the same color 
class, thus reducing its available color space by only one 
color. 

The algorithm to determine which vertex 1S next, 
considers the rows in the adjacency matrix for all the pre= 
viously colored vertices, and adds the columns. The aighest 
column total then indicates which vertex is to be colored 
next. Previously colored vertices can be eliminated by sub- 
tracting an arbitrarily high number from the running column 
total, for each vertex, when it 1s chosen. In the case of 
ties, the vertex of highest degree is chosen. 

S. Least colors available (LCA) 

In all the previous static orderingsSs an assumption 
was made that the best vertex to color next, at any point in 
the trees, was independent of the specific color asignments 
made up to that point. We were simoly trying to find the 
vertex that was ‘most orobably’' a better choice than the 
others. In this dynamic orderings when expanding each node, 
we actually choose the vertex which will result in the 
fewest children being created, thus minimizing the branching 
factor. Since a child node is created for each available 
color which the next vertex can be colored, we must consider 


all uncolored vertices in the path to the node being expand= 
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ed and determine which vertex has the fewest colors avail= 
able. To do this note that ideally we would like to choose 
a vertex which is adjacent to at least one vertex in each 
color class used so far. Thus a new class would have to be 
created and there would be only one child node created. 
Short of this ideal situations we would like to choose the 
vertex which is adjacent to at least one vertex in most of 
the color classes so far created. | 

This 1s accomplished by creating a new matrix CC 
similar to the adjacency matrix such that there is a row in 
CC for each color class generated and a column for each ver 


tex. Then we define 


cc... = ! if (vr-j) © E for some vertex v in color class 


iS 0 otherwise. 


i iff vertex j cannot be colored color i. 


ie © eur og 


The matrix CC is easily created from the adjacency 
matrix A as follows: Let ce. be row i in CCy and a. be row 3 
] 7 
in A, both represented as bit strings. 


initialize al} cc. . = 0 

for i from 1 to NUM=COLOR=CLASSES 

for each vertex j in color class i 
cc. = ce. V a. 


The next step is to generate a vector called SUM such that 
cc. 1f j 18 uncolored 
ae 1 J 


-1 if j is colored 


1 
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Then the value of j for which SUM is a2 maximum is 
the vertex which has the fewest color classes available to 
which it can be added. Again, for consistency, ties are 


broken by choosing the vertex of highest degree. 


E. CONTROLLING ALGORITHM 

Inasmuch as the purpose of this research was to experi 
ment with various orderings of the verticess the flow of 
control was fairly straightforward. After only a few trial 
cases it was discovered that the primary limitation to the 
size of the graph which could be colored was computer 
storage spacer, not time. As the graphs got larger, more 
nodes were required in the search tree until the available 
memory was used up at about 4000 nodes. Although unsigned 
variables were used in the node structures and space was 
conserved wherever possible, this limitation could not be 
Significantly changed. For even the best orderings, this 
meant graphs of about 30 vertices or less. For empirical 
data then, graohs were generated with from 10 to 30 vertices 
and edges from 20% to 80% of the maximum possible. For each 
Qraph size, 100 graphs were generated (some of which were 
provably identical and many which were at least isomorphic) 
and each ordering method was used to color the graph. The 
graph was then divided at its articulation points and each 


biconnected subgraph was colored using each ordering funce 
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tion. Data about the size and shape of the search tree was 
collected and tabulated (for all the above Samples 
(appendix 8). 

Note that 100 trials of each graph size is a relatively 
smal} sample, but was a necessary restriction to keep the 
entire program execution time to within a few hours. 

The following algorithm shows the setup used for obtain= 


ing the experimental] data. 


algorithm MAIN 


initialize 
For n from minsize to maxsize (n is number of vertices) 
begin 
compute maxedges = (n) (nel) /2 
For fullness from 30% to 70% of maxedges 
begin 
compute number of edges e = fullness * maxedges 
For trials from 1 to 100 
begin 
generate random graph (n vertices, e edges) 
For each ordering method 
begin 
color the graph 
record data (expansions, nodesr branching factor) 
divide the graoh 
For each biconnected component of the graph 
begin 
color the component 
record component data 
end 
sum the component data 
end 
end 
compute averages of all data collected 
print out data 
end 
END 
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The function MAIN in appendix A is the controlling algo= 
rithm for the entire process. 

F. GENERAL COMMENTS ON ALGORITHM DEVELOPMENT AND THE 

PROGRAMMING LANGUAGE '‘C! 

In all the programming done during this research, ef 
forts were made to save time and storage space whenever 
feasible. The first attempt to design a coloring algorithm 
used a backtracking method which 18S common for tree 
searches. However, the size of the search tree was so large 
that the storage required to stack the variables used in 
each recursive call to the backtracking function was exorbi-= 
tant. For that reasons an iterative scheme was chosen and 
implimented. 

Further savings were made by using dit strings and il bit 
variables whenever possible. Since the maximum graph size 
was to be less than 3e vertices, the long integer variables 
of 32 bits in 'C' were ideally suited for the adjacency mao 
trix and many other needed variables. Also 'C' has a good 
assortment of bit manipulation operators including multiple 
shiftSs offs ands and exclusive or. Using these led to time 
SAa8VINQS iN many cases. 

Space was conserved by using the unsigned variable type 
im afrayS and records. With this convention, the length of 


the storage cells reserved were only as long as needed. For 


54 





example, the vertex numbers from 0 to 31 could be stored in 
only 5 bits. Another oracticer not necessarily considered 
good programming technique was the fairly extensive use of 
global variables. Though risky, it did save the time re- 
quired to pass parameters and the space to store local varis 
ables. 

Overall, this author feels that 'C' is a very appropris 


ate language for research of this type. 
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ives ResulL iS 


A. THE PARAMETERS 

In attempting to determine the relative ‘goodness’ of 
several computer programSr one must first find some measur= 
able parameters to be minimized or maximized to indicate the 
best orogram. For approximation algorithmsr one 1s very 
often not Interested In how long a program takes as long as 
it is some linear function of the problem size or possibly 
even O(nlign). What one measures is the closeness of the apo 
proximation to the actual optimal solution. 

For this research, the solution to the graph coloring 
problem was to be the exact chromatic number tn every case. 
The coloring algorithm was to find the chromatic number of a 
given graph and the ‘goodness’ of the vertex ordering 
methods was to be determined by how fast the algorithm ran 
and how much space was used. The specific questions to be 
answered were: 1) Is there a best*time ordering for all 
Qraphs?, 2) Is there a best-space ordering for all graphs?, 
3) Are they the same?, 4) What tyne, or sizer, of graphs 
favor which ordering if there is a difference?, and 5) Is 
one ordering better on the average, but Snotiber better for 
worstrcase situations. 

The results of the many exoeriments which relate to 


these questions is tabulated in appendix 8. The discussion 
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here is of a general nature and will summarize the results. 

There were basically two separate categories of measurements 

made = those related to time and those related to space. 

Though they are closely related in this problem, the parame= 

ters measured are distinct and will be discussed separately. 
1. Time 

In determining the difference in time=complexity of 
the 5 coloring methods, there are 2 modules of the program 
to consider = ORDERNODE and COLORALL(Gof which EXPAND is a 
part). These functions are the only two whose times= 
complexity 18 a function of the vertex ordering as well as 
the graph size. 

For the 4 static orderings the time=complexity is 
easy to determine. For consistency, each ordering merely 
filled an array called next, where next(i) was the vertex to 
follow 1 in the ordering. Random is therefore clearly O(n) 
since the vertices are ordered by vertex number. One 
traversa| through the array next is all that is needed to 
fil? it. DEGREE is 0(n*) since the adjacency matrix has to 
be traversed to find the degree of each vertex, then a sim 
ole sorting of the vertices is O(nIign) so the overall come 
plexity is 0(ne). DVEC 3 is O(n) if the number of itera- 
tions is kept constant(such as 3 in this case). It would be 


conceivable to iterate n times for the best possible order= 
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ing in which case DVEC mn would be 0(n°). ADJACENCY is also 
Bin’) because a running total was kept of the relative voosi= 
tion of the unselected vertices as each new one was deter-= 
minede Thus the adjacency matrix was addressed in only one 
row for each new vertex. For the dynamic ordering LCA, no 
work was done in the function QRDERNODE but rather in the 
function PICKNODE which was called every time a node was to 
be expanded. PICKNODE itself has a time complexity of n3 


operations to create the n x n matrix of 


since it takes n 
color classes CC and each element in CC requires O(n) to 
compute. But the bigger concern is that the time=complexity 
of the branching funetion EXPAND is n3 times as great for 
LCA than for the static orderings. 

For the static orderings the time-complexity of gen 
erating the search tree has 2 components = the time to ex= 
pand a node and the number of expansions needed to find a 
solution. For LCA, a third components, the time to find the 
next vertex, must also te considered. 

Once the next vertex is founds, the time required by 
the function EXPAND is O(n). The only iterative part of 
that function is following the path to the root and then 
Generating a child node for each possible coloringe Since 
the depth of the tree can be up to n=l and the number of 


children can be ny, this function is O(n). 
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So far, all the time complexities have been polyno- 
mial. If this were aiso true of the number of expansions, 
then one could claim that the overall algorithm would run in 
polynomialetime Since the product of any number of polynomio 
als is also a polynomial. Unfortunately, that is not the 
case for worstecase situations. Consider for a moment the 
size of the full search tree generated by this algorithm if 
no pruning 1s allowed. Refer to figure 11 which shows the 
first four levels of the search tree. The notation xry in 
each node of the tree indicates that the node represents an 
assiaqnment of color y to vertex x. Note that at the root 
there is 1 node _ followed by 2 nodes at ltevel 2. Then at 
level three there are 5 nodes, 2 children of one in level e 
and 3 for the other. To determine the number of nodes at 
level k, recall that one color class is added for each level 
in the tree. AlSoOr a node at level k represents one of the 
ways of partitioning the vertices into k or fewer color 
classes. Thereforer the number of nodes at level k is the 
number of ways to partition a set of n elements into k= or 
fewer classes. 

The Stirling numbers S(nrm) reoresent the number 


of ways to partition n elements into om distinct 
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classes({Ref. 12) and are generated by the following recure 
rance relations: 


Stn tliykJe=s ome kel) + k * S(n-k) 


i 
for all on 
l 


To understand this relations, consider the following: 


where S(n,n) 


and S(n,1) 


Given that on elements can be partitioned into k classes in 
S(nek) wayS», there are two ways to add the (ntl)st element. 
Fiesty, it can be added to any one of the k classes for each 
way of partitioning n elements. The number of ways of doting 
this is given by k * S€ nek). Secondlysr, the n elements can 
be partitioned into k=l classes and then the (nt+l)st element 
can be a single element in a new class. This 1S the reason 
for the S(n-ke1) term. 

But the number of nodes at level k is the sum of the 
Stirling numbers S(nem) for m from 1 to ke. These numbers 
are called the Bell numbers 8 after E.7T.Bell (Ref. 13). B 
1$ the number of partitions of a set with n elements into 
any number of distinct classes up to ne That is 


EP = > S€n¢k) 
K 


To see how fast these numbers grow, the first 8 are: 
1,225715,52,203,877,4140 
In other words, at level 8 in the search tree there 


would be 4140 nodes and the total number of nodes in the 
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tree would be the sum of B. for 1 from 1 to 8. To see what 
the upper and lower bounds on BY are, consider again the 
tree in fig. 11. The fewest number of children of each node 
is found in the leftmost path down the tree. In this path, 
every node has two children. In the rightmost path down the 
tree, a node at level +? has i children and this 1s the larg=- 
est number of children for any node at that level. Thus the 
total number of nodes at any level i lies somewhere between 


the two caseS»s, or 


Since BY is growing faster than ZA it 1S easy to see 
that as n gets larger the possibdble size of the tree 18 going 
to far outweigh any of the polynomials so far discussed. 
Ooviously, unless a significant amount of oruning is done, 
the time to expand every node at each increasing level mn 
the tree is going to grow far too rapidly. Since the search 
Strategy was fixed for all orderingse the most significant 
way to compare the time savings of each ordering was to 
count the number of nodes which needed to pe expanded. 
Therefore, the number of calls to the function EXPAND was 
used as the primary measure of the time complexity of the 


various orderings. 
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a SD ace 

There are only two data structures in these programs 
whose size is dependent upon the vertex ordering: the array 
of nodes in the search tree and the array of pointers used 
for the priority queue. But the number of nodes in the 
Queue iS always less than the total number of nodes inthe 
tree since expanded nodes are removed from the queue. Since 
the number of expansions hasS been recorded for ae time 
analysiSs only the total number of nodes in the search tree 
at the time an optimal solution is found need be recorded 
for a space analysis. 

In order to see the shape of the search treer anoth= 
er parameter, the average branching factors was recorded for 
each level in the tree. The average branching factor is de= 
fined to be the average number of children of each node 1n 
the search tree. A full binary tree would thus have. an 
average branching factor of 2.0. If the branching factor 1s 
recorded for each level in the search trees, one can see the 
‘shape’ or the ‘'narrowness' of the tree. The most desirable 
Situation iS the one in which the branching factor Stays at, 
or as close as possible, to 1.90 in the first Several levels 
of the tree. This indicates that the maximum pruning took 
Dlace aS early in the search as possible, thus reducing the 


search soace by more nodes for each prune. The "'goodness' 
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of the ordering functions can thus be seen by how long the 
branching factor stays at 1.9 and how low it stays when it 


does start to increase, 


B. TIME ANALYSIS 

In analysing the time required by the program for each 
ordering method there were two cases of interest: average 
time, and the time required in the worst case. Figure i2 is 
a plot of the average number of expansions required for each 
ordering method when graphs containing 30% of the maximum 
edges were colored. 

For the worst case graphs the natural logarithm of the 
number of exoansions was plotted since exponential growth 
was expected. Figure 13 is a plot of the worst case graphs 
for 30% edges. These plots are not smooth and in fact are 
quite ragged because each point actually represents” the 
coloring of one particular graph in the sample of 100. Note 
that the ordering methods parallel each other in the dips 
and hills of the plots indicating that the same graph was 
usually the most difficult to color for all of the methods. 

The significance of these two plots is in the total dom- 
inance of the LEAST COLOR AVAILABLE method of ordering the 
vertices. Fhough all the other methods eventually ‘explod= 
ed' even in the average case, LCA stayed very close to 


linear for the full range of graphs up to 3e vertices. 
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Furthermore, the worst case analysis shows the number of ex- 
pansions for LCA leveling off between 20 and 30 vertices. 
It would be an interesting experiment to carry the graph 
sizes out past 32 to see if and when LCA also starts to show 
explosive growth. 

Another results though not Shown 1Nn graphical form 1s 
that the time required by any of the ordering methods 1s 
lowest for graphs with very few edges as wel] as for graphs 
which are close to being complete. For any given value of 
ne empirical evidence indicated that the most time 1s rem 
quired when the number of edges is in the vicinity of 50% 
(see appendix B). This is an intuitively obvious result if 
one considers what happens in the search tree.w With very 
few edges, very few colors will be needed so the number of 
children of a node will be minimal. The Searches, in this 
case, will tend to proceed down the left side of the search 
tree. Most of the right side will be pruned away because 
the cost of the nodes in the right side of the tree will be 
greater than the cost of the optimal solution. 

In the case where there arfe many edgesr a new color will 
be required at each level and the search will tend to stay 
very close to the right side of the tree. The left side 
will be pruned away since many of the nodes on that side 


will represent illegal colorings due to edges. 


o7 





mee oOPACE ANALYSIS 

Just as for time, the space analysis is of both average 
and worst case situations. Figure 14 is a plot of the aver- 
age numver of nodes in the search tree when the first op- 
timal solution was found. Figure 15 is a plot of the natur- 
al log of the maximum nodes needed for any of the i100 
graphs. Again, the number of edges is 30% of the maximum 
possible. 

Not surprisinglys the space analysis results are very 
similar to the results found in the time analysis. In every 
case LCA is superior in the sense that it uses less space 
and its growth rate is much improved over the other methods. 

Since the relationship between the number of expansions 
and the total number of nodes is the branching factor, it 
comes as no surprise that LCA exhibits the lowest branching 
factor at the top of the search tree. Figure 16 shows the 
average branching factor as a function of depth for each of 
the 4 heuristic orderings. Note that LCA and ADJACENCY have 
the same branching factor for 3-4 levels but then ADJACENCY 
grows faster in the midrange. Note also that both of these 
methods have a higher branching factor in the middle of the 
tree than either DEGREE or ODVEC 3. This fact shows how 


coowerful 3 or 4 prunings at the top of the tree can be. 
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D. ADVANTAGE OF DIVIDING THE GRAPHS 

In the first stage of this coloring process, the graphs 
were colored just as generated with no attempt to divide 
them into biconnected components. In the second stage, they 
were divided and then each component was colored and the 
number of expansions required was Summed up for all the com- 
ponents. For sparse graphs where there were many very smal] 
components, the results of this experiment were quite sub- 
stantial. However, aS the graphs became more dense, the 
original graph was biconnected and no Savings were derived 
from dividing. 

Figure 17 shows the size of the largest biconnected com- 
ponent as a function of the number of edges. In this plot, 
the ratio of the size of the largest component to the origi-= 
mal graph size is plotted for 4 cases. Note that since the 
maximum number of edges grows as né , when the number of 
vertices in the graph increaseSr the graph becomes bicon= 
mected with a smaller percentage of edaqes. Thus it iS rea 
sonable that the empirical evidence gained here indicates 
that a graph with 100 vertices is biconnected 99% of the 
time with 6% edges whereas a graph with 50 vertices is 
biconnected 99% of the time when it has 12% of the maximum 


vertices included. 
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In investigating the effects of dividing on graph color= 
inar note from fig. 17 that with 30% edges, the largest 
biconnected component is very close to the original graph 
for all ne Thereforer in order for dividing to be benefi- 
cial, random graohs must be very sparse. For the empirical 
results shown in fig. !18 graphs were generated with from 10% 
to 30% edges. fhe benefit of dividing is then shown by 
plotting the ratio of the number of expansions needed for 
the divided graph to the number of expansions for the undio= 
vided graph against the percentage of edges oresent. 
Inasmuch as the relative benefit was the same for ail order- 
ing methodss only LCA was used for this graph. 

It appears as though the benefits from dividing are ino- 
Significant since only very sparse random graphs have arti 
culation points. But many reale-life poroblems which are 
modeled by graphs may be quite dense yet have articulation 
points. This is true particularly in transportation = and 
communications problems where some of the vertices may 
represent junctions or hubs through which al] traffic must 
flow. Consider the city of Chicago in a nationwide tran- 
sportation network for example. In these cases graph divid= 
ing may result in a much faster coloring. 

Also there are poroblems which have graph models”) which 


are indeed very sparse. Consider, for example, the exam 
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scheduling problem in a large university. The vertices of 
the grapn are the classes and there is a edge between two 
vertices if at least one student is enrolled in both 
classes. The question then 1s: How many exam periods are 
needed to schedule all final examinations with no conflicts? 
Consider the hypothetical but realistic case of a university 
with 200 classes, and a average of 20 students per classe 
Since many students take the same curriculum together it is 
not inconceivable to have a graph with an average degree of 
each vertex of 5S. This would result in a graph of 200 ver= 
tices and 500 edges or 2 1l/e %. This iS a very sparse graph 
which is likely to have several articulation points. In 
fact, the biconnected components are likely to be the vari- 
ous departments of the university. Dividing the graph prior 
to solving this problem would surely result in a time’ and 


space savings. 
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IV. CONCLUSIONS 


The branch and bound algorithm described in this paper 
can be used to find the chromatic number of graphs. Though 
no empirical evidence was gained for graphs larger that 32 
vertices, it is felt that somewhat larger graphs can also be 
colored in a reasonable amount of time on the average with 
the LEAST COLOR AVAILABLE method of ordering the vertices. 
Further, some improvement in the lower bound function or 
search strategy might be possible which would make the algo 
rithm more efficient. For example, some lookahead procedure 
might conceivably improve the method of assigning costs to 
each node of the search tree. Alsor the priority of a node 
in the queue was primarily the number of colors in the par= 
tial coloring and only secondarily a function of the depth 
in the tree. A possible improvement is to make the priority 
some algebraic combination of the two factors thus giving 
more emphasis on the depth in the tree. This might result 
in some faster colorings. Since the emphasis here was oan 
finding the best branching function, little effort was ex- 
pended in looking for improvement in the other 2 components 
of the branch and bound strategy. 

For the five branching functions investigated, LCA seems 
far superior for all graphs in both the time required and 


the amount of computer storage space needed. Though al] 
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five functions have worst case time and space complexity 
which apoears to be exponentials the LCA method has an avers 
age case time complexity which appears to be nearly linear 
in the range of graphs from 5 to 3e vertices. 

Though the research and experiments done here focused on 
the exact algorithm for graph coloring, a few very minor 
changes would result in an approximation algorithm. The 
function EXPAND could be changed such that it only generates 
the first possible child node rather than all children of a 
node. No priority queue would de needed and the number of 
nodes required in the search tree would just be ny, the 
number of vertices in the graph. So the soace-complexity of 
the algorithm would be linear in new The time required would 


be O(n? 


) if LCA were used for an ordering functions but it 
would be interesting to see how good its ability to approxi- 
mate the real chromatic number 1s compared to other order= 
NGS. 

Dividing a graph at its articulation points prior to 
coloring is an effective way to reduce the time required to 
find the chromatic number. However, unless the graph is 
known to have articulation points or unless it is very 
sparse, the effort spent in finding the biconnected com- 


eponents may not be worth it. When a graph does have articus 


lation points, there is a tworfold benefit gained by divid- 
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ing the graoh. The first 187 of course, that 
the resulting graph 18 smaller and thus easier 
The second benefit comes from the fact that the 


components tend to be dense which makes coloring 


19 


the size of 
to, color. 
biconnected 


easier. 


variables. 


This appendix contains the source 


puter programs used 


The remainder of the 
grouped together according to their 


The modules and the function 


MODULE 


EXTERNE DECS 


MAIN 


GENGRAPH 


ORDERNODES 


COLORALL 


DIVIDE 


TREEMANAGER 


in this research. 


APPENDIX A 


into various functions which Share some 
These variables are listed 


modules contain 


FUNCTIONS 


none 
main() 
gengraf() 
random() 
ordernode() 
picknode() 
Go1orala() 
expand() 
convo() 
divide() 
Groc() 
getnode() 
addq() 


gettop() 
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Phestina stor 


The program 


in the module EXTERN 


in this appendix are: 





Jfxrnrekw kw &® ke ® wk eR wR &® wR wR wR wR &e Re Re RR RR Rl lll 
GLOBAL MACRO DEFINITIONS 

meRek ke we wR we me le lll le le RR mR ke me ke Re ke me ke OK ke lk kf 

#define NUMNODES 4000 /x Maximum number of tree nodes */ 

#define NUMQ 2500 /x Maximum size of the priority Q x*/ 

define NULL 0 /x The empty pointer k/ 

define NUMTRIALS 100 /x Tre number of trials for each size 

of graph x/ 
tdefine MAXSIZE 30 4x The largest graph size (n) x/ 
#define MINSIZE 10 4x The smallest graph *x/ 


J/xnnwrenewerke ke we we &e we eR we me eR lw le lll lll Kl Kl Kl le Re KR 
EXTERNAL DECLARATIONS 
mek we we we we we Re me me Re me me me mR Ok kkk kk ke KS 
/x trnode 1s the structure of the nodes 
which make up the search tree */ 
Struct trnode { 


unsigned depth :8; 
unsigned cost 3:8; 
unsigned vertex 38; 


unsigned color 3:8, 
struct trnode *parptrye 
d? 
/x tnode 1s a array of search tree nodes. 
Gis an array of oointers to tree nodes */ 
struct trnode tnode (NUMNODES] - xq f[NUMQ) ; 
/xk ndata 18 a record of data regarding the 
graph vertices to be used tn ordering the 
vertices by the branching function. 
order is an array of records */ 
Struct ndata { 


unsigned next 57 
unsigned degree "57 
unsigned flag > oi; 
} order(3e); 
int fd; 4x file descriptor for output */ 
int nz 4x number of vertices */ 
int count? /x number of calls to EXPAND x/ 
int hinode; /x first vertex to color. */ 
int index; /x index of next available node in the 


array of tree nodes = also the number 
of nodes in the search tree at any time *x/ 


int qsize,;z /x index into the priority queue */ 

int chrom, /x The chromatic number of the graph *«/ 
int color{32); /*x The color class of each vertex */ 
int ordnum; /x The ordering function being used */ 
int Saven; 4x The saved value of n for DIVIDE */ 
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int grafl32) (32) /*x The adjacency matrix */ 

float numsubs; /x*x The number of biconnected subgraphs */ 

float bigsudn;s /* The size of the largest biconn. como. */ 
float bigsube /x Number of edges in bigsubn */ 

long int graph{32]); /*x The adjacency matrix A in bit strings */ 
int sonct(3e)3 /*x A counter of child nodes for each level */ 
int parct({32e); /*x Parent counter for each level */ 


i a a a ee a ee ee ee a a a a ee ee a a a a a a a, a ee | 
MAIN 

wh we wk Ow we Ow Om Ow OR OU OK Oe OK OR Ok Ok Ok OK Ok kk ke ke ok ke ke ke Oe kK 
/x This function provides for the flow of 
control throughout the program and records 
the data which 1s to be collected. It also 
contains the code for any output generated */ 

main () { 


int full; /x The percentage of edges */ 
int edges, /x The number of edges */ 
int isjes /x jteration variables */ 


int ndflagl6)s /*x Too large of a search tree flag for each 
ordering method */ 

int trial; /x The trial number for each graph size */ 

float maxedges; /*x The maximum edges possible for graph */ 

float maxobrl6j; /*x The maximum expansions required */ 

float maxnd(l6}); /*x The max nodes required */ 

float avgorl6); /* The average expansions required */ 

fioat avgndl6); /* The average nodes in the search tree */ 

char *label (6); 


4x The output file and labels */ 
fd = copen("/work/rerauten/out 3",'a');3 
label {1) = "RANDOM"; 


label (2) = "DEGREE"; 
label (3) = "DEG 3"; 
label (4) = "ADJCEN"; 
label (SJ = "LCA"; 


/x. For each value of n desired */ 
for (n=MINSIZEsn<=MAXSIZEsnt+=5) |¢ 


/* Compute maximum edges possible */ 
maxedges = (n * (nel))/2;3 


/*x For each value of the percentage of 


edges desired */ 
for(full=eszfull<=207 fullt+=2) ¢{ 
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/x compute number of edges */ 
edges = full * maxedges/100; 


/x Print headings */ 
printf (fd,"0 =23d, EDGES =%dd(%2d%%)0, 

Nn,edges, full), 
orintf(td,"%5s%10s%10s0,"order","“branches","nodes"); 
printtlfds"K415sziles"-"avqg max","avg max "); 
printf (fd,"0); 


/x Initialize data collection variables */ 
for(islzi<sSzitte) { 

maxbr [ji] 0; 
maxnd (ij 
avgor (i) 
avgnd{i] 
ndflag(i)] = 0, 
) 


6 tb 
= 


/x For the number of trials desired *«/ 
for(trial=l;trial<=NUMTRIALS;trialt+t+) { 


/x Generate a random graph */ 
gengraf(n,edges); 


/x For each ordering method */ 
for(Cordnum=S;ordnum<=Szrordnumt++) { 
if(! ndflaglordnum)) { 


/x Divide the graph and color the components */ 
if(! divide()) ndflagflordnum] = 1; 


4x Record data */ 
avagbrlordnum) += count, 
avgondlordnum] += index; 
if€count > maxbe lordnum) ) 
maxberflordnum) = counts 
iffindex > maxndflordnum) ) 
maxndlordnum) = index; 
) 
} 
} 
/x Print out the data collected */ 
for (izlsi<=S;,i+t) { 
if(! ndflagfi]) { 


/x Compute the averages */ 


avgorli] /= NUMTRIALS; 
avgndli] /= NUMTRIALS; 
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pervect (ta, «68S",labehli}); 
printft( fd, *Z4.0f25.0f ",avgbrli] smaxbrli)); 
orintf(fde-"%60.0f%5.0f "“pavandlil ~maxndli) ); 
Oernttcrar, 0), 
} 

else printf(fd,"%6s NOT ENOUGH SPACE AVAILO, 

label fil); 
} 


} 
cexit(); 


} 


i a a a ee a a ee a ee ee a a a a a a ee Se ee i 
GENERATE GRAPH 
meek ke we ke we he me Re UR le lm OR OK Om km Ok OR OH ek ke ke ke KS 
/x This function generates a random graph with 
nm vertices and e edges and records it in a 
global adjacency matrix */ 


gengraf(n,e) int neers { 


INt irvjrarbes 
long int convo (); 
/x initialize the adjacency matrix */ 

for (1=0;7i1<=n7itt+) { 

graphli] = 0; 

for(js=0; j<e=n7 j++) graflid Cj) = 0; 

grafl(il (i) = i; 

} 
while (e > Q) { 

/x Get e random vertices */ 
random(n); 
random(n); 

/x Find the next ungenerated edge */ 
while(graffa}) (b) == 1) { 

a=(aZn) + 13 

if (a == 1) b = (6b A mn) + I, 

) 


a 
b 


/rx Record the new edge */ 
graffal(b) = 1; 
graffb) (a) = 1; 
==, 
) 
/x Put adj matrix in array of bit strings */ 
for (i=lsi<=enz,ite+) §{ 
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graflil fli) = 0; 
for(jztsj<an-jtt) iff€graflil (j)]) graphli) t= convolj); 
} 


ak &£ & ke RK KK kK eR Ke Ker ehlUu hlUhlU hUCU LhlUCU hl KhlUlU hl lhl hl hl hl hl! hl RR 
RANDOM 
a a ee ee a ee ee a a a a a a, A 
f/x This function is used to get random numbers 
between 1! and the input argument on. 
taken from Grogono [ref. 2] */ 
int random(n) int nz { 


static long int seed 18403; 
int 2 


seed = (25173 * seed + 13849) % 65536; 
i= seed 4 Nn; 

return(itli); 

} 


Jxenw kk ke &w & ww wR we ® wR we wR ke wk we wh we we wR eR ke we me le lm le le le 
ORDERNODE 

ewe k we we we we me lm le le Om OK Oke OK Ok kk kkk lk kek ke ke ke ke ks 
/x This function orders the vertices of 
the graph for each of the static ordering 
methods. It does so by filling the array 
order.next with the vertex to follow each 
vertex in the graph. */ 

ordernode() { 


emt irsjrks /k jteration variables */ 

int save; /x the saved last vertex */ 

int temphi-, /x the temporary next vertex */ 

int max? /x any maximum needed to maintain x*/ 


int colsum(32); /* used in Adjacency ordering to keep 
track of number of oreviously colored 
vertices a vertex is adjacent to. */ 
int maxcol; 4x The column with the max 1's in colsum */ 
iNt hideg,s /x The highest degree of any vertex */ 
float dvec(4) (32); /*x The vector of degrees in DVEC 3 x/ 
float maxdvec; /x The max vector calculated in DVEC 3 x*/ 
long int masks /* A bit manipulation mask */ 
long int convo (): /* A converting function (SUPP MOD) */ 


max = QO; 
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f/x initialize the vertex data */ 
for(i = OF i<=enz7i1¢++) ¢ 
order(i]).degree = 0; 
orderlil.next = 0; 
orderli]l flag = Q; 
mask = grapnli); 


/x Compute the degree of each vertex */ 
for(j=lsj<ens jer) { 
if (€(mask & 01) == 1) orderli]).degree += 1; 
mask >>= 1; 


} 


/x find the max degree and the vertex */ 
ifforderli).degree > max) { 
max = orderli) .degree; 
hinmode = i- 
} 


4x Order the vertices depending on the 
current ordering number */ 
Switch (Cordnum) { 
case 0: 
break; 


Jx RANDOM x/ 
case 1: 
for (i=0z7i<nzitt) orderli) next = i + 17 
hinode = 1; 
break; 


/x by DEGREE */ 
case 2:3 
order(0]).next = hinode; 


/«x .flag indicates already colored vertex */ 
order lhinode) .flag = 1; 
Save = hinode; 


/x For each vertex find the next one */ 
forl(izlri<snzite) { 
max = Q; 


4x Find the highest degree of the uncolored 
vertices */ 
for(j=Hlrj<enz jor) { 
if€(Corderlj] -degree > max)&&(! orderlj).flag)) { 
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max = orderlj).degree; 
temphi = 373 
) 

) 


/x record the next vertex */ 
ordertlsave] .next = temphi; 
orderltemohi].flag = i; 

Save = temphi; 
) 
break; 


Je by DOVEC 3 x/ 
case 3:3 


/x initialize Oth iteration to degree */ 
for(iztlsi<znzit+) dvec(0) (i) = orgerli) -degree; 


/x For 3 iterations */ 
for(islszi<=s3ss7it¢te+) f 


/*x For each vertex */ 
for(jaltsj<snz jee) { 
dvecl(ijJ (j] = 0; 


4x For every other vertex in the graph */ 
for(ksi;zk<snjzk+rt+) { 
mask = convo(k);. 


/x If they are connected compute summation 
for the ith iteration */ 
if(€(graph(j)] & mask) > 0) 
dvecfiJ (j] += dvecli-1) (k); 
} 
} 


4x Order the vertices */ 
4x The ordering method 18 identical to case 2 
exceot for the ordering parameter */ 
Save = 0; 
forlislri<senzi+t+) { 
maxdvec = 0; 
for(j=alisj<enz jee) { 
if((dvec(3) (j] > maxdvec) &8& ($ orderljJ].flag)) { 
maxdvec = dvec(3) (jj; 
temohi = j?3 
) 
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order(save] .next = tempohi; 
order(temphil].flag = 13 
Save = tempohi, 
} 

hinode = order([0) .next; 
break; 


4x by ADJACENCY x/ 
case 43 
/x initialize x*/ 
for(islszi<snzitt) colsum(i) = 0; 
order (0) .next = hinode; 
order lhinode).flag = 1; 
save = hinode; 
colsum(save} = 33; 


4x For each vertex find the next one */ 
for(jalsj<enelsje+) { 

max = Q; 
mask = O1; 

/zk tncrement the bit string colsum with each 

vertex added to the set of colored vertices */ 
for(li=szlsri<=enz7zit++) { 

if((graph(save) & mask) > 0) 
colsum{i] ++, 


iIf€colsumli] > max)max = colsum(i); 
mask <<= 1; 
} 


/x Order the vertices = again the same 
ordering method 1s used as 1N case @ except 
the ordering parameter is by value of colsum */ 
hideg = 0; 
for(i=zltsri<enszi+t) { 
if(Cecolsum(li) ==max)&& (order (i) .degree>hideg) 
&&(! orderli).flag)) { 
hideg = orderli).degree; 
temphi = 17 
} 
} 
order [save] .next = t 
orderl(temphil.flag = 
colsum[ftemohil] -= 33 
save = temphi; 


break; 


88 





famoyeteEAot COLORS AVAILABLE *x/ 


case 53 
/x no work done in this function. 
See PICKNODE x/ 
break? 


Jxnrm tk &® & wk we & & ww ke wR Re RH w we eR wR we we Re OR Ue Rl RR Ke Ue 
PICKNODE 

x *« ®* k ke eR we ke RH ke Re RelURlUKRlUlUlURlU Rll Rll lll Rl 
This function is called from the routime 
EXPAND to get the next vertex to color. 
If one of the static orderings 1s being 
used it just looks up the next vertex in 
the table order.next. If the dynamic 
ordering LCA is being used, it must trace 
the search tree to the root to determine 
the next vertex to color. x / 


picknode(nodeptr) struct trnode xnodeptrs { 


/x Declarations x/ 
long int ccl32e) -vdonersliong); 
int maxrisumek,i,suml(3e) ,nodenum-shideg; 
struct trnode *tempptr; 


f/x If static ordering, just table lookup */ 
iffordnum < 5S) returnlorder (nodeptr=>vertex} next)? 


f/x initialize 
vdone 1s a bit string to record which 
vertices have already been colored. 
cc 18 the color class matrix 
sum counts the colors not avail for each 
vertex, */ 
vdone = 0; 
forl(izlri< 


=n7z itt) { 
ec{i] = Q; 
sum{(i] =0; 
) 


temooptr = nodeotr; 


/x Trace from node N to root. */ 
while (temoptr != NULL) { 


/zx fi111 color class matrix */ 
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ccltempptr=>color) t= graoh(tempotr= vertex] ; 


/x record colored vertices */ 
vdone != convbo(temootre=>vertex); 
temootr = tempeotr=>oparptr, 


/x For each of the n vertices */ 
for(k=lek<enzk+er) { 


/zx If the vertex k is not colored */ 
if((vdone & convb(k)) == 0) { 
isum = Q; 


/*x count the color classes to which vertex 
kK cannot be added */ 
forlisztszi<s=nodeptr=>costsitt) 
1SUM += 


{ 
long] & Ceclij>>Ckel)); 
} 
sum(k] = 


1SuUM; 


/x keep track of fewest number of color 
classes available */ 
if€isum > max) max 


=> 1SuUM; 


Je Of all 


the vertices with the fewest color 
classes available, find the one with the 
highest degree */ 
for(isltri<=enzi+e) { 
if((sum(i}) == 


== max) && (i€convo(i) & vdone))) { 
ifCorder{i]) .degree > nideg) { 

hideg = orderli).degree; 
modenum = 137 


} 


/x Return the vertex number of the next 
vertex to color. */ 
return(nodenum); 

} 
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SS ee ee ee i i i a a a a a a a a 


CULOR-ACE 


——x~ xk kX XH uw Ee KH HR EH KH KR KE RR KRH RR RR KR RR ReCUlRCUCURCR RH RR 


This function creates the root node of the 

Search tree, puts it in the priority queue 

then calls the function EXPAND to expand 

the top node in the queue until the top of 

the queue iS a complete coloring. */ 
colorali() { 


struct trnode *nodeptr,*gettoo(); 


int 2 
/x Initialize for new coloring */ 
k = 07 
index = Q; 
qQqsize = 0; 
count = Q; 


tnodel0).cost = 313 
for (€12=031<NUMQ;i++) gli] = &tnodel(0); 


/x Get an empty node, fill it with the 

data for the root of the tree and 

out it In the queue. */ 
if €! (nodeptr = getnode())) return(0); 
nodeptr=>parptr = NULL; 
nodeptr=>deoth = 1; 
nodeptr=>cost = 1; 
nodeptr=>vertex = 
nodeptr=>color = 1 
addq(nodeptr),; 


hinode; 
, 


/x ‘forever! x/ 
for(s7) { 
/x Get the too node in the queue */ 
nodeptr = gettop(); 


/x Tf this is complete coloring, quit */ 
if(nodeptr=>depth == n) break; 
else if (! expand(nodeptr)) return(0); 
} 


4x Record the chromatic number */ 
chrom = nodeptr=>cost; 


/x Record the color class assignment of 
each vertex in this path. */ 
while (nodeptr=>parptr != NULL) { 
color (nodeptr=>vertex] = nodeptre=>color; 
nodeptr = nodeptr=>parptr,; 


a1 





} 
color (nodeptr=>vertex] = nodeptr=>color, 
return(1); 
} 


I i 
EXPAND 
x *x*ikk Re ke & & & & kw & & & & wR RH ke HR & RH RHR wR RR we Re we we ® 
This function expands the node pointed to by 
the variable nodeptr. It creates a child 
node for each possible color available as 
well as a node for one new color. */ 


expand (nodeptr) struct trnode *nodeptr; { 


struct trnode *tempotr, *newnode, *getnode(); 
mnt vertexnumei; 
long int mask,colavail,convo(); 


/x increment counters */ 
count *#= 1; 
parct [nodeptr=>depopth]) ++; 


/x Get the number of the next vertex to color */ 
vertexnum = picknode(nodeptr); 


/x create new nodes for all colorings using 
only existing colors */ 
colavail) = O1; 
for (i=l7i<nodeptr=>costritt+) { 
colavail] <<= 1; 


colavail i= 01; 


/x colavail is a bit string with a one in the 
position of every available color. */ 
} 
tempoptr = nodeptr,; 


/x For each node N in the path to the root */ 
for (373) { 


/x Tf Nevertex and vertexnum are connected */ 
if((graphi[vertexnum) & convol(tempptr=>vertex)) $= 0) { 
mask = convb(temoptr=>color); 


/x Remove N.ecolor from colavai!) */ 
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if (Ccolavail & mask) !2= 0) colavail f= mask; 
} 
if (temoptr=>parotr == NULL) break; 
temopotr = tempptr=>parptr; 
} 
/x Tf any colors are left */ 
if (colavail != 0) { 
mask = O11; 


/x For each color available x/ 
for (i=l; i<=nodeptr=>costritt+) §{ 
if (€C{colavail & mask) != 0) { 


/x Generate child node *x/ 
if ((€newnode = getnode()) == 0) return(0); 
newnode=>depth = nodeptr=>deoth + 1; 
newnode=>cost = nodeptr=>cost, 
newnode=-vertex = vertexnum, 
newnode=>color = 1, 
newnode=>parptr = nodeotr,; 

4x yncrement child counter */ 
sonct [nodeptr=>depth) ++; 


/x Put the node in the queue */ 
addq(newnode); 
) 
mask =<< 1; 
b 
} 

/x Generate one more node with a new color */ 
if ((newnode = getnode()) == 0) return(0); 
newnode=>-dJepth = nodeptr=r-depth + 1; 
newnode=>cost = nodeptr=>cost + 1; 
newnode=>vertex = vertexnum,s 
newnode=>color = nodeptrercost + i; 
newnodes>parptr = nodeotr,s 
sonct (nodeotr=>depth] ++; 
addq(newnode); 
return(1); 

} 


Jxnnwnee hm we we we we we mR mR Ome me lm Uwe lm ll Oe lO ke OOK 


SUNVER TT? 10 (61 } 


keke kK wk Re ew Oe Ke OR we ke Oe Oe Oe KR OR OK Ok ke OR OK OR kk ke ke ek 
/x This function converts an integer argument 
between 0 and 16 to a bit string which has a 
1 in the position of the argument and O's 


93 





everywhere else x*/ 
long int convob(arg) int args { 


int 2 
long int bitpos?s 


bitpos = O13 

forl(izltri<argritt+) bitpos =<< 1; 
return(bitpos); 

} 


Jxrxerek kw &® & & & & & we we R K wR wR & wR wR eR we wR wR we eR eR eR KR RR 
DIVIDE 

wk ke we eR we em eR Re Re RR RR OR Owe ke Ke he OR Ok kek kk kk tke oS 
/x This function divides the graph into its 
biconnected components. I[t then creates a 
new adjacency matrix far each component and 
calis the function COLORALL to color the 
component. It records the data for each 
component and sums it up for all comps of a 
graph. When it returnse the sum data is in 
the variables set up by MAIN, x*/ 


divide() { 

int number(3e); /zx to record the visiting order of verticesx/ 

mat back {32}; /x record the lowest vertex to which a back 

edge exist */ 

iNt stackv(3el]; /x A stack of the vertices */ 

mt num, /x To maintain the next number x/ 

Int wrv,too, /x vertices = top is top of stack */ 

int stvptr; fx pointer to top of vertex stack */ 

int steptr; /x pointer to top of edge stack */ 

int nume; 4x number of edges 1n biwcomponent */ 

int numv; /x mumber of vertices in dizcomponent */ 

mt vnum[(3e); /x array to record the vertex numbers as 
the edges of the edge stack are popped */ 

int totcount; /x count the total calls to EXPAND for al) 
components */ 

int totindex; /x count the total nodes used */ 

int table(3e); /* to convert the bicomponent to a graph */ 

Int tableptr, /x pointer into table */ 


int vertli,verte; /x the vertices of the edges in new graph */ 
long int convb();3 
int i,j} 
/x An edge in a biconnected component defined 
by the two vertices at its end */ 
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struct edge { 
unsigned vl: 8; 
unsigned ve: 8; 
} stacke[350); 
/x Stacke is an stack of edges */ 


4x Save off the oringginal value of n = the 
biconnected components will be different */ 
saven = nz 


4x Initialize al) variables x*/ 
number(l] = 1; 
num = 23 
stackv{0}] = 1; 
top = ls 
for(v=erv<=ensvet) numberlv] = 0; 


steptr = ll, 

stvetr = 0; 

totcount = 03 

totindex = 0; 
/x Begin a depth first search of the graph 
exploring forward edges */ 

forwards: 


/x While there is an unexplored forward edge 
from the top vertex on the stack to some 
vertex w */ 
while((w = proc(top)) != 0) { 
/x Stack the edge (toprw) on the edge stack */ 
stacke(++steptr}]).vl = top; 
Stacke(steptr]).ve = wi 
4x If w has already been explored give the min 
of number(w) and back(top) to back(top). */ 
if€number(w]l > 0) { 
if€numberlw) < back{top]) back{top] = numver (wl); 
} 
/x Else stack w and explore it */ 
else { 
number(w)] = numtt; 
back{w] = numver (wl); 
Stackv{(++stvptr] = w, 
top = we, 
} 


4x If there i8 no forward edge then backup 
to the last vertex which stil! has a forwara 
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edge which is unexplored */ 
backups 
if(stvptr > 0) { 
f/x try the next to top vertex */ 
v = stackv(stvptr - 1)7 
/x if it has a number less that the back of 
the too vertex then we have an articulation 
point and the bi connected component 18s on 
top of the edge stack */ 
iff€back[ltop) >= numberl(v]) { 

/x increment and initialize counters */ 

numsubs ¢= 1; 

nume = 0; 

numv = QO; 

for(i=isi<=saven;sit+tt) graphli)] = 0; 

tableptr = 0; 

table(0) = 0; 

for(i=0;7;1<=savensitt+) vnumli] = 0; 
4x pop edges until end of component */ 

while((number(stackel(sSteptr) .vi}] >= numoerlv)) && 

(number (stackel{steptr].v2]) >= numper{v} ) 
&& (steptr >= 0)) { 


Jz increment data counters */ 
numertt; 
if(++vnum(stackelSteptr].vijJ = 
if(t++vnum(stackel(Steotr).v2] = 
verti = 0; 
verte = 0; 

/x check to see if either vertex is already 

in the table = if sor, get the index which 

1S going to be the new vertex numoer */ 
for(iz=07i<=tableptreit+t) { 


= 1) 
= 1) 


Numvest 7 
NumMvete 7 


if(tableli] == stacke(steptr]).vl) verti = if 
if(tableli] == stackxelsteptr].ve) verte = i; 
} 


/zx if not in table, put it there and get new 
vertex number *x/ 
if(verti == 0) { 
tableptrtt; 
table(tableptr) = stackel(steptr) .vl; 
vertl = tableptr? 
) 
if(verte == 0) { 
tableptrtrt; 
table({tableptr] = stackel[steptr] .ve; 
verte = tableptr; 
} 
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/x record the edge between the new vertex 
numbers */ 
graohi(vertil) } 
graph{vert2) ; 
/x decreme 
steptr=e; 
} 


= convb(vert2d); 
= convol(verti); 
mt edge stack pointer and do again */ 


/x The biconnected component iS now recorded 
in the adjacency matrix and ready to color */ 
/x Is this the biggest component of this graph */ 
if(numvy > bigsubn) bigsubn = numv; 
if(nume > bigsube) bigsube = nume; 
/x new nis size of the table */ 
n = tableptr, 
/x Call the ordering function */ 
ordernode(); 
4x Color the comoonent */ 
if€! colorall()) return(0); 
/x Increment the cummulative counters */ 
totcount += count, 
totindex += index; 
} 
/x Tf no arte pte pop the vertex stack and 
explore any forward edges of the vertex 
which 1S Now ON too */ 
else if(back{top] < back{v]) back{v] = backlitop); 
top = stackvl=-stvotr] ; 
goto forward; 


} 


/x When al! the connected vertices have 
been explored, check for another unconnected 
component and explore it */ 
for(i=27i<=saven;sitt) { 
if(number li] == 0) { 
number li) = 1; 
num = 2; 
stackv{0) = i; 
top = i, 
stvptr = 0; 
Ssteptr = =-l; 
goto forward, 
} 


/x When everything 18 done and all components 
are found and colored, set the counters to the 
cummulative totals, reset nmr and reset the 
adjacency matrix */ 
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eount = totcount, 
index = totindex; 
Nn = savens 
for (i=lLFi<=enzitte+) { 
graffi) Ci] = 0; 
for (jslrjy<sn7 jet) 
if (graf(i) (j) 
if (graffi) (j] 
} 


“ues 


2) Geattii ti] = 1; 
1) graph(i}) := convo(j); 


} 
return(1); 
} 


fxknere ke ® ® ® ® wR wR wR wR wR wR wR we R wR we we me me Ke OR RR KR UR 
PROC 

kk we we ke Oe OR me me OR OK OR OK Om Om OR Ok OR Ok Ok OR Ok ke ke ee ke ke kK 
/x This function determines if we can 
oroceed in the depth first search exploring 
forward edges. It marks the adjacency 
matrix with ae as each edge 1s explored 
so that no edge is explored twice */ 

proc(i) int is { 


int j7 


/x Is there a vertex adjacent to vertex i 
On an edge which hasn't been explored. If so 
mark it explored and return the vertex number */ 
for(j=17j<=saven; j++) 
if(graflij fj} == 1) { 


graffi) (j) = 2; 
graff{jl (iJ = 2; 


return(j); 
} 
return(0); 
} 


a a a ee 2 ee 2 ee 2 
TREE MANAGER 
ze we wR we me wm OR Kl lm lm ORO he Om Rk OR hk ke Om mk kk OK ke oS 
4x This module contains the functions 
necessary to manage the search tree nodes 
as well as the priority queue. */ 
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fx wre ke & wk we we kK kK kR we & KR RK KR kK Ke Ke KR Ke eR Re eR wR Re Re RK 
GETNODE 

ee 2 a 2 a 2 2 a 2A 
/x This functions maintains the array of 
nodes and returns a pointer to the next 
available node when called. Index is a 
global variable x/ 

struct trnode xgetnode() { 
/x jindex gets initialized in COLORALL *x/ 

index += 1; 

if Cindex < NUMNODES) return(&tnodelindex] ); 

else return(0); 

} 


f/x xk ke ke ke RR ke &® we we we me well KhlUhlu hl hlhKhlh hl hl Khel lhl KOK Oh 
ADD TO QUEUE 

ewe kkk wk ke we Kw Ke Ke KR Khu he KR OK UK OK OKO kkk ke kek KS 
/x This functions adds a new: node to the 
priority queve. It puts it in the last 
position and then ‘sifts it up’ until it is 
in its proper location */ 

addq (ptr) struct trnode *ptr; { 


struct trnode x*temp; 
int tej? 
/x Put the new node at the bottom */ 

qsize += 1; 
if (qsize == NUMQ) printf(fd,"Q OVERFLOWO),; 
alaqsize)] = ptr; 

/ex j 18 the node being sifted up */ 
asize; 


} 
4x 1 18 the parent of j */ 
m= j/2;3 
/x Nhile j is not the root */ 
while(j > 1) { 
/x compare priority of j to i */ 
iff{alil=->cost < qlj]*>cost) return; 


else if (qlile>cost == qlj)=>cost) { 

if CalilJ=>deopth >= alfjl=>depth) return; 

} 

/x Exchange nodes and continue */ 

temp = ql{j); 
qalj] = qlil; 
qalijJ = temp; 
j) = i; 
i = j/e; 
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fix ® &* & &R ke eR eK KR RR RR RR ell Kell RRR RR RR KR 


GET TOP 


mk x wk ke kK eK KR OK ClURClUKhlUhlUhlURhlU hl hl hlURhlU hh hh KR KR Ke OK OR OR KS 


/x This functions returns a pointer to the 
tree node which 1s on top of the priority 
queue. The queue 1S a heap. When the root 
of the full binary tree iS removed, the node 
in the last position 1S moved to the root and 
then ‘sifted down’ */ 

struct trnode *gettop() { 


Struct trnode *saveptr,*temp, 
int ivje, 
/x save off the root node */ 
saveotr = qlilj; 
/x move the last node to the root */ 
all} = qlasizel),; 
/x decrement the aqsize */ 
qsize == |i; 
/x sift the root node down */ 
/x 4 18 the parent node being looked at */ 


/x j is the left child of i */ 


/x While were not out of the heap */ 
while(j <= qsize) { 
if (j < aqsize) { 


/x Get the child of 1 with least cost */ 
if (qljlercost > ql{jtlle>cost) j = jl; 


/x if both the same costes get greater depth */ 
else if (qljl=>cost == q{jtl] ->cost) 
if€(qljle>depth < qlj+tler>deoth) jf = jet? 
} 
/x Tf both children have lower priority, quit */ 
if (qlil=>cost < q{jle>cost) break; 
else if (alile>cost == qljle>cost) 
if (qlile>deopth >= qlj]*>depth) break; 


/x exchange node i with jlits highest 
Priority child) x/ 


temp = q(jl; 
qa(jjJ = qlil; 
gli) = temo; 


/*x get i = to j and do it again */ 


=! 
ib 
em 
=e 
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/x when done, return the save off root */ 


ipn(saveptr); 
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APPENDIX 8B 


This appendix contains output listings from the various 
programs run to generate the plots described in Section IV. 
The first listing is the data for graph coloring without dio 
viding for graphs with from 10 to 30 vertices. The edge 
values are 30%, 50%, and 70%. 

The second Jisting 18 the output from the DIVIDE routine 
which divides the graohs and records data regarding the size 
of the biconnected components and the number of components. 
The graphs ranged in size from 6 to 100 vertices. 

The next listing contains data from tne routine to 
divide the graohs and color the components. Graohs with 25 
vertices and 2% to 14% edges were used in this experiment. 

The last listing shows the effects of the amount. of 
edges when nis held constant. For this experiment n = 20 


while the edges vary from 20% to 80% by 5% intervals. 
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ne ee a ee a, a ee a, a a a, a, a ee, ee, ee a a, a, ee, ee, a, ee, a a a, a, 


COMPARISON OF 5 ORDERING METHOOS FOR SEVERAL GRAPH STZES 


x & & &® &e & eR we Re eRe ee we Re ee Rell lll lel lll OK 


m= i0, EDGES = 13(30%) 
ORDER BRANCHES NODES 
avg max avg max 
RANDOM 19 65 se 6115 
DEGREE 10 14 21 29 
OVEC 3 9 1i el en 
ADJCEN 9 9 el 24 
LCA 9 9 21 24 
BRANCH FACTOR BY LEVEL 
1 2 3 4 5 6 if 8 a 10 
mammom 1.55 1.92 1.96 2.14 1.75 2.11 1.85 2.43 2.01 0.00 
Meeree 1.25 1.59 1.82 1.94 2.29 2.41 2.57 2.78 2.88 0.00 
Myeco 5 1.04 1.45 1.96 e@.22e 2.45 2.47 2.70 2.83 2.67 0.900 
ADJCEN 1.00 1.29 2.16 2.44 2.41 2.45 2.68 2.83 2.89 0.00 
men 1.00 1.29 2.15 2.2c9 2.47 2.61 2.1/4 2.87 2.92 0.00 
N= 10, EDGES = 22(50%) 
ORDER BRANCHES NODES 
avg max avg max 
RANDOM 29 i1e9 53 239 
DEGREE 10 18 el 35 
OVEC 3 9 ol 20 25 
ADJCEN 9 11 20 eo 
LCA 9 9 20 25 
BRANCH FACTOR BY LEVEL 
1 2 3 4 5 6 7 8 9 10 
meer 1.54 1,81 1.72 1.99 1.63 2.24 1.57 2.13 1.75 0.090 
Meer 1.08 1.36 1.59 1.76 1.93 2.25 2.57 2.66 3.13 0.00 
Mees 1.01 1.18 1.42 1.97 2.27 2.42 2.56 2.73 3.13 0.00 
SeyeeN 1.00 1.00 1.35 2.13 2.44 2.50 2.65 2.80 3.16 0.00 
mee. .00 1.00 1.35 2.07 2.35 2.45 2.7e 2.87 3.24 0.00 
N= 10, EDGES = 31(70%) 
ORDER BRANCHES NODES 
avg max avgG max 
RANDOM 33 160 56 250 
DEGREE 10 17 19 eT 
DVEC 3 9 16 18 26 
ADJCEN 9 10 18 ee 
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BCA 9 9 18 ee 
BRANCH 
1 2 . “ 

RANDOM 1.08 1.50 1.48 1.98 
DEGREE 1.01 1.12 1.24 1.36 
mee 5 1.02 1.05 1.19 1.21 
ADJCEN 1.00 1.00 1.00 1.13 
men t.00 1.00 1.00 1.13 


N= 15, EDGES = 31(30%) 


BRANCHES 
avg max 


NODES 
avg max 


ORDER 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


oe 
17 
LS 
14 


48 
38 
36 
55 


88 
4e 
34 
ee 


173 
8e 
66 
44 


FACTOR BY 


5 
hay) 
ieoe 
eo) 
leer 
1.78 


NOT ENOUGH SPACE AVATIL 


BRANCH FACTOR BY 


l 2 5 4 

RANDOM NOT CALCULATED 
Meeeee loci 1.59 1.67 1.86 
meee 5 1.06 1.24 1.77 2.07 
ADJICEN 1.00 1.02 1.64 2.13 
men 1.00 1.02 1.62 2.07 


li le 13 14 

RANDOM NOT CALCULATED 
DEGREE 2.66 2.83 3.05 3.44 
OVEC 3 2.77 2.88 3.20 3.43 
BevcEN 2.78 2.98 3.15 3.44 
mee 2-82 3.07 3.32 3.60 


N= 15, EDGES = 52(50%) 


ORDER BRANCHES 


avg max 


NODES 
avG max 


RANDOM 
DEGREE 
BYEC 3 
ADJCEN 

LCA 


24 
LS 
16 
1D 


72 
47 
61 
33 


47 
38 
55 
a2 


Leo 
89 
115 
66 


15 


0.00 
0,00 
0,00 
0.00 


NOT ENOUGH SPACE AVAIL 


eee 


eau? 
1.98 
2.15 
24 
Bees 


DEVEL 


2.03 
2.28 
2.48 
2.4) 


0.00 
0.00 
0.00 
0.00 
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7 
1.43 
Co 
a. 51 
2629 
Sele 


2e09 
ero! 
2.47 
ce. 44 


0.00 
0.00 
0.00 
0.00 


8 
Coes) 
2.65 
2.73 
Cet 
2.76 


eace 
AD 
2.46 
Ceol 


0.00 
0.00 
0.00 
0.00 


9 
1.48 
3.10 
5.14 
5.14 
Bs 


2.56 
2525 
2.49 
ao? 


0.00 
0.00 
0.00 
0.090 


10 

0.00 
0.00 
0.00 
0.00 
0.00 


10 


2.49 
2.43 
2.65 
eel 


0.00 
0.00 
0.00 
0.00 





RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
MyEC 3 
ADJCEN 

LCA 


N= 


ORDER 


RANDOM 
DEGREE 
BVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
BVEC 3 
ADJCEN 

LCA 


a = 


ORDER 


15, 


20, 


BRANCH FACTOR BY 


1 2 3 4 
NOT CALCULATED 
fee f50 1.47 12.57 
$1.04 1.17 1.34 1.59 
1.00 1.00 1,04 1.49 
1.00 1.00 1.04 1.47 


le i 14 
NOT CALCULATED 
eesce ceo ec. / i 5,14 
Pit eed ec. /0 S516 
BOE sO5 Te Ot o.ec 
Pou es!0 2.93 35.38 


11 


EOGES = 73(70%) 
BRANCHES NODES 
avg max avQ max 


5 


lie 
ine 0 
eave 
7 


{> 


0.00 
0.00 
0.00 
0.00 


NOT ENOUGH SPACE AVAIL 


75 
56 
se 
28 


hes 
96 
66 
35) 


ab 
38 
31 
29 


BRANCH FACTOR BY 


! 2 3 4 
NOT CALCULATED 
me te 0! Jets 1.351 
Weve t.0c 1.06 1.24 
1.00 1.00 1.00 1.00 
m0 1.00 1.00 1.00 

i Ve 13 14 
NOT CALCULATED 
eas ¢«na e€,/0 3.18 
eee! 2646 2,68 3.15 
2.45 2.54 2.85 3.14 
eet? 2-695 2.91 3.44 


EDGES = 57(30%) 
BRANCHES NODES 
avQ max avG max 


5 


1.4e 
1.56 
1.08 
1.08 


1 


0.00 
0.00 
0.00 
0.00 


BEVEL 


1.81 
1.96 
2.43 
e.cl 


0.00 
0.00 
0.00 
0.00 


LEVEL 


o0 
1.60 
1.49 
Lea? 


0.00 
0.00 
0.00 
0.00 


A, 


Vee) 
2.04 
Cro 
eee l 


0,00 
0.00 
0.00 
0.00 


tei 
Ll.// 
2.03 
leon 


0.00 
0.00 
0.00 
0.00 


eau 
eet 
2.41 
esa? 


0.00 
0.00 
0.00 
0.00 


ied 
2.01 
Coon 
2.14 


0.00 
0.00 
0.00 
0.00 


2.30 
ewer 
eae2 
2eot 


0.00 
0.00 
0.00 
0.00 


Leos 
2.05 
2.40 
2.14 


0.00 
0.00 
0.00 
0.00 


10 


2.25 
2.45 
2.41 
Cie 0 


0.00 
0.00 
0.00 
0.00 


10 


aecv 
2555 
ea 
Cea5 


0.00 
0.00 
0.00 
0.00 





RANDOM 
DEGREE 
myEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
OVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


N = 


ORDER 


RANDOM 
DEGREE 
PVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
mcC 
ADJCEN 

LCA 


20, 


NOT 
63 
33 
eo 
ee 


601 
Ceo 
418 
247 


SNe 
118 
208 
lel 


te 
70 
37 
$3 


ENOUGH SPACE AVAIL 


BRANCH FACTOR BY 


1 e 3 a 
NOT CALCULATED 
Poss? ere 1. Se 
Peco el.o4+ 1,08 1.87 
1.00 1.00 1.49 1.84 
1.00 1.00 1.49 1.77 

11 Me 13 14 
NOT CALCULATED 
2.17 2.30 2.41 2.49 
eeomee a! 2. 58ue.5) 
Bel 0 26.94 €.58 2.50 
Bee 2-04 2650 2.56 


EDGES = 95(50%) 
BRANCHES NODES 
avg max avG max 


5 


Vig tele. 
2.01 
ewe 
2.07 


15 


ewok 
aya 
2 ol 
EoD 


NOT ENOUGH SPACE AVAIL 


110 1345 203 2592 
76 1206 143 2308 
46 881 95, leo 
31 648 Sou heod 


BRANCH FACTOR BY 


1 2 3 4 
NOT CALCULATED 
melo 1.54 125/ 1.55 
Peo? 1.26 1.56 1.60 
1.900 1.00 1,00 1.14 
1.00 1.00 1.00 1.14 

11 ee 13 14 
NOT CALCULATED 
eo? e¢.0e 2.91 2.38 
Peloee.co 2€.15.e.e7 
eo ee DS 2.50 2.20 
Ee? 2.56 2.44 2.33 


5 


leo 
1.66 
loo 
or 


15 


2.41 
2.45 
2.40 
2.42 


EEVeEL 


6 


1.60 
2.05 
2uee 
ewe 


16 


2.56 
2.68 
2.57 
ears 


EEVeL 
6 


9S 
een 
2.00 
1253 


16 


2.54 
2.45 
222 
ee 


106 


3 


1.76 
2.06 
2.42 
2.56 


Pea ible 
Cet 
2.81 
Ceo 


Ze 


lef 
1.63 
2il 
1.9e 


17 


2.66 
2.66 
2.61 
2e0f 


2.08 
2.28 
2.40 
2c4> 


18 


2.80 
a7! 
2.86 
5.18 


1.59 
io? 
ae 
2-15 


18 


Cad 
Ce 
eo 
5.09 


2.08 
2.41 
ero 0 
ee ot 


19 


Deo 
3.26 
3.31 
3.50 


1.64 
1.93 
2.25 
e.c/ 


5.04 
5.05 
3.11 
oe) 


10 


eell 
2.32 
a. 1G 
Csoe 


eV 


0.00 
0.00 
0.00 
0.00 


tet) 
1.90 
Cee 
Ge ae 


20 


0.00 
0.00 
0.00 
0.00 





N = 20, EDGES = 133(70%) 


ORDER BRANCHES NODES 
avg max avQ max 


RANDOM NOT ENOUGH SPACE AVAIL 
DEGREE 52 18e Se oe? 
mec 3 36 «119 Sve ecud 
ADJCEN et 4} 43 76 

mea 19 35 41 66 


BRANCH FACTOR BY LEVEL 
l 2 3 4 5 6 7 8 7 10 
RANDOM NOT CALCULATED 
Meenee 1,02 1.03 1.07 1.17 1.45 1.45 1.67 1.64 1.88 1.87 
mee > 6f.00 «12601 «1.05 1.07 1.25 1.56 1.57 1.75 1.82 1.52 
ADJCEN 1.00 1.00 1.00 1.00 1.00 1.01 1.13 1.63 1.97 2.30 
Meta c0 12.00 1,00 1.00 1.00 1.01 t.te 1.63 1.88 2.04 


RANDOM NOT CALCULATED 

Meee 1.90 1.67 2.17 2.29 2.52 2.64 2.81 3.10 3.32 0.00 

OVEC 3 1.98 2.11 2.29 2.43 2.58 2.72 2.84 3.07 3.54 0.00 

ADJCEN 2.42 2.52 2.64 2.80 2.65 2.76 2.99 3.13 3.47 0.00 
Meee ec) Co4) €,535 2.65 2.77 2.98 3.19 3.37 3.91 0.00 


N = 25, EDGES = 90(30%) 


ORDER BRANCHES NODES 
avg max avqQ max 


RANDOM NOT ENOUGH SPACE AVAIL 
DEGREE NOT ENOUGH SPACE AVAIL 
Mec 5 126 1158 ealweccesdo 
ADJCEN 63 411 tsSce 019 

men 35 115 eee 


BRANCH FACTOR BY LEVEL 
1 2 3 “ 5 6 7 8 9 10 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
fees 1-20 £.35 1.58 1.85 1.80 1.93 1.98 2.05 1.79 1.84 
Meee 1.00 1.00 1.32 2.06 2.12 2.05 2.20 2.09 2.15 2.07 
eer. 00 1.00 i.32 1.97 1.91 2.08 1.98 2.03 2.08 2.08 


11 ae 13 14 15 lo 17 18 19 20 
RANDOM NOT CALCULATED 
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DEGREE 


NOT CALCULATED 


mee 5 1.7/7 1.85 2.00 2.97 2.95 
MmceEN 2.09 2.03 2.09 2.15 2.26 
men 2.02 2.035 2.99 2.ee 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


N = 


ORDER 


RANDOM 
DEGREE 
DVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
pyEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
OVEC 3 
ADJCEN 

LCA 


RANDOM 
DEGREE 
OVEC 3 
ADJCEN 

LCA 


25, 


22 25 24 
NOT CALCULATED 
NOT CALCULATED 
ead C€.70 2.05m5.47 
pee O10 9 6 II te SO 
2c? 5.08 3257 3.85 


21 


EOGES = 150(50%) 
BRANCHES NODES 
avg max avg max 


Cn) 
25 
C0 6 


0.00 
0.00 


NOT ENOUGH SPACE AVATIL 
NOT ENOUGH SPACE AVAIL 


178 1328 340 2568 
97 932 193 1829 
SiGe) 59 ome i] 


BRANCH FACTOR BY 


1 e 3 - 
NOT CALCULATED 
NOT CALCULATED 

Me09 71.28 1.29 1,49 
1.00 1.00 1,00 1.12 
1.00 1.00 1.00 1.12 


11 te 13 14 
NOT CALCULATED 
NOT CALCULATED 
ioe iece 1,87 1.80 
meer 1698 2.902 1.92 


eove2s01 1.95 1.95 


ei 22 23 24 
NOT CALCULATED 
NOT CALCULATED 
pe 6.9? 2,06 5.05 
2.535 2.368 2.59 3.08 


Ges? se eOs 2.94 35.29 


> 


lao 
1.40 
1.34 


is) 


ico 
2.04 
ieee 


25 


0.00 
0.00 
0,00 


Cleo mest 663) 2651 
eeOme cree IO eC. 5c 2.5 
2.42 2.48 2.54 2.61 2.75 


LEVEL 


etd, 
1.93 
1.06 


16 


1.9¢ 
lee 
2.04 
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les 5 
2.18 
et 


wi 


1.83 
2.26 
2.08 


1.87 
2.24 
ono 


19 


20 


2.16 
ec 
eat 





Reve 
5 6 7 8 9 10 


Mer eS, EDGES = 210(702%) 
ORDER BRANCHES NODES 
avg max avg max 
RANDOM NOT ENOUGH SPACE AVAIL 
DEGREE NOT ENOUGH SPACE AVAIL 
OVEC 3 NOT ENOUGH SPACE AVAIL 
ADJCEN 87 723 169 1431 
LCA 34 122 64 225 
BRANCH FACTOR BY 
1 e 3 4 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
DVEC 3 NOT CALCULATED 


meocEN 1,00 1.00 1,00 1.00 
LCA 1.00 1.00 1.00 1.00 


11 ine PS 14 


RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
OVEC 3 NOT CALCULATED 


MOUCEN eC.ece 2.54 2.22 2.16 
mee i.78 2.05 2.t2 1.95 


el ee 23 24 
NOT CALCULATED 
BEGREE NOT CALCULATED 
EywEC 3 NOT CALCULATED 
MemeeN 2.10 2.44 2.59 2.85 
mea 2.56 2.58 2.84 3.41 


RANDOM 


N = 30, EDGES = 130(30%) 


ORDER BRANCHES 


avg max 


NODES 
avG max 


RANDOM 
DEGREE 198 
OVEC 3 88 412 eS 
ADJCEN 63 691 134 1391 

LCA 31 57 iy 


1504 


1.00 


120 


1.48 
1.43 


ie 


Ort Otol ait > 


iS Lo 17 18 19 20 


ZelS 
wiomewWeges0oue.Ud 2.10 2. 


NOT ENOUGH SPACE AVAIL 
Seoreded 
815 


BRANCH FACTOR BY LEVEL 


1 2 3 4 


5 6 7 8 a 10 
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NOT CALCULATED 

le? ee Oe 
1.64 1.74 
ies. oO 
ode 1.65 


RANDOM 
MeGREE 1.36 1.73 
mveC 3 1.21 1.48 
ADJCEN 1.00 1.00 

LCA 1.00 1.00 


le 13 14 
RANDOM NOT CALCULATED 
MeeReEe 2.09 2.18 2.08 1.97 
mec 5 1.95 2.05 2.16 2.33 
feeeN 2.28 2.21 2.25 2.20 
mea ce iD 2.14 2.23 2.19 


11 


el ee 23 24 

RANDOM NOT CALCULATED 
MeoeREE 2.26 2.59 2.435 2.42 
Mee > Cw 41 2.32 2.54 2.53 
MeceN 2.41 2.42 2.42 2.58 
mee 2246 2.91 2.92 2.60 


N= 30, EDGES = 217(50%) 
ORDER BRANCHES NODES 
avg max avg max 


Lo? 
1.94 
Ce 10 
1.94 


15 


Cau) 
eae 
een 
e.cl 


25 


Gat 
eae 
Ciena 
ele 


RANDOM NOT ENOUGH SPACE AVAIL 


DEGREE NOT ENOUGH SPACE AVAIL 
PEC 5 289 1825 Soe >olc 
ADJCEN 79 1995 163 3959 
BEA 32e 76 70 148 
BRANCH FACTOR BY 
1 2 3 “ 2 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
mee > 61.08 1.ie 1.33 1.59 1.70 
MeeereN $.00 1.00 1.00 1.01 1.13 
mee 1.00 1.00 1.00 1.01 i.,i2 
Li lie 13 14 i> 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
Meee > 61.82 1.80 1.76 1.94 1,89 
ADJCEN 2.30 2.29 2.33 2.27 2.28 


eer 
let 
Cecl 
2.04 


lo 


2.04 
ewer 
2.3e 
2ecl 


26 


2.42 
2.54 
515) 0, 
ome 


EEVEl 
6 


er? 4 
1.64 
1.63 


lee 
1.92 
eeco 
en 0 


Nd 


ele 
2ec4 
e.cl 
2.c4 


eT 


2a 
ee 
ae) 
S60 5 


1.89 
ce 01 
1.88 


17 


1.86 
vo. 
Geee 
Aa 


18 


CS 
2.356 
2.16 
@.5¢ 


ae 


on 
Ciao 
eno 
3.12 


ealS 
2.20 
2.00 


18 


2.18 Celt 
Caro 2.e8 


es 
ie 7c 
e2oc 
2.13 


19 


eas 
Gel 0 
a 
2.56 


29 


2.80 
2.84 
2.87 
5.41 


leo? 
22> 
e201 


19 


1 9S 
e05 
2.25 
eto 


20 


ee Os 
2.34 
Coad 
2.38 


30 


0.00 
0.00 
0.00 
St, 


10 


beoe 
ae) 
2.08 





men 2618 2.09 2.15 2.20 


ae ao 24 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
OveC 3 2.33 2.43 2.44 2.51 
ADJCEN 2.40 2.37 2,43 2.45 
men 2.50 2.55 2.66 2.64 


en 


N = 30, EDGES = 304(702%) 


NODES 
avg max 


BRANCHES 
avG max 


ORDER 


RANDOM 
DEGREE 
Mee 3 215 1519 
ADJCEN 46 1003 

LCA 30 67 


394 2896 
96 2003 
65) ice 


2.15 
eo 
2ea4 


ao 
2.80 


NOT ENOUGH SPACE AVAIL 
NOT ENOUGH SPACE AVAIL 


BRANCH FACTOR BY 


1 2 3 4 

RANDOM NOT CALCULATED 

DEGREE NOT CALCULATED 
Myec 5 1.00 1.00 1.04 1.05 
ADJCEN 1.00 1.00 1.00 1.00 
LCA 1.00 1.00 1.00 1.00 


ie iS 14 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
meme > 1.97 2.03 2,04 1.99 
mesceN 1.6/7 1.97 2,19 2.35 
mee 1.66 1.88 1.96 2.09 


11 


el ec 23 24 
RANDOM NOT CALCULATED 
DEGREE NOT CALCULATED 
eee 5 2.35 2.49 2.64 2.70 
BeyeeN 2.60 2.69 2.89 2.81 


mem c.68 2.70 2.86 2.98 


5 


1.08 
1.00 
1.00 


> 
io5 
2.48 
2.19 
eS 
2.84 


2.93 
Rare) ae 


2.18 
26 
2.62 


2.0e 
2.9e 


LEVEL 


yee 
1.00 
1.00 


lo 
1.72 
eS PY 
2.26 
26 
2.98 


3.04 
ee 


Cie 7 


a) 


2.835 
22/8 


3.06 


1.36 
1.00 
1.900 


17 
1.68 
ene 
2.34 
ee, 
5615 


3.18 
3.48 


2.24 


28 


2.94 
2276 
So 


1.34 
1.00 
1.00 


18 
Zaut 
2.45 
2.54 
28 
3.34 


Bro? 
>! 


5.17 
Dre 
Sit 


le? 
ae 
ac 


Hh) 
Bie 
eo 
2.45 
eo 
Sea 


3.83 
4.23 
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OUTPUT FROM GRAPH DIVIDING ROUTINE 


zk &e ek eR eR eR Re melRelURlUelUKRlUell lll le le lk 


0.00 
0.00 
0.00 


10 


1.65 
lnc 
1.28 


ey 
eo 
2.55 
2.45 
30 
0.00 


0.00 
0.00 
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AVG SIZE LARGEST SUB 


AVG # SUBS 
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N 


edges 
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18.0 


i 
13.8 
13.6 
15.9 
14.0 
149.0 
14.0 


4.1 
10.7 
lS<c 
15.8 
15.9 
i>. 7 
16.0 
16.0 


4.9 
10.5 
toe? 
15.9 
16.9 
vec 
17.6 
aE 
17.9 
ies? 
18.0 
18.90 
18.0 
18.0 


4.6 
era 
16.¢ 
ie.) 
19.1 
19.6 
19.8 
19.8 
20.0 
20.0 
20.0 
20.0 
enn 9 
20.0 


28 ac 
piles 5 
34.6 
37.9 
41.0 
44.0 
47.0 


Gel 
ia.8 
22.8 
Geel 
51.9 
36.9 
40.0 
44.0 


4.9 
13.4 
20.9 
Clee 
32.7 
S67 
41.3 
45.5 
49.9 
Dioeo 
DO. 0 
62.9 
66.0 
609.9 


4.6 
16.0 
24.7 
De 0 
>e8.0 
44,3 
a7 
594.8 
a0 
65.0 
67.79 
75.0 
80.0 
85.0 
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fee? 
16.5 
18.9 
202 
ei.l 
21.4 
Cit 
21.8 
21.8 
eles? 
Gere 0 
e2.9 
e2.9 
e2.0 


5.2 
lSe¢ 
19.7 
ec.) 
aaa 
Cael 
23.8 
Con? 
24.0 
24.0 
24.0 
24.0 


eee! 
io. 


bo. 
24.7 
32.6 
40.0 
45.5 
50.6 
S10 <) 1 
61.5 
66.4 
Pleed 
77.0 
82.0 
87.0 
92.0 


Sec 
19.7 
50a 
39.3 
46.4 
53.4 
59. 
65.6 
72.0 
78.0 
B42.0 
90.0 


See 
tome 
30.0 
38.5 
46.5 
Deo 
60.6 
Sle 
73.0 
79.7 
85.9 
91.6 
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eae > 
26.6 
e7.e 
27.6 
orien 
le? 
27.9 
28.0 
28.0 
28.0 
28.0 


56 
17.6 
259 
Coe 
28.1 
28.7 
29.5 
255 
eta! 
C969 
C959 
30.0 
30.0 
30.0 
30.0 
30.0 


6.2 
31.7 
37.8 
39.4 
39.9 
40.0 
40.0 
40.0 


Ole 
38.7 
46.7 
48.5 
49.5 
49.9 
50.0 
50.0 
50.0 


15.05 
Sa. i 
o1.3 
69.1 
76.3 
aoc 
90.9 
98.0 
104.8 
112.0 
119.0 


5.6 
Geal 
56.0 
44.8 
Bie 
SSG 
7o. 5 
78.e 
S355 
92.8 
99.7 

107.0 
114.0 
Bega. 0 
128.0 
135.0 


6.2 
49.0 
(9 
98.6 

119.8 
140.0 
160.0 
180.0 


665 
61.2 
7D. oO 

120.8 
148.7 
174.9 
200.0 
225.0 
250.0 





70 


80 


90 


100 


60 


120 
re 
180 
210 
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70 
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140 
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280 
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350 
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120 
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240 
280 
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90 
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180 
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360 
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450 
495 
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6.7 
46.1 
56.0 
58./ 
59.7 
59.9 
60.0 
60.0 
60.0 
60.0 


6./ 
a5 7 
64.6 
68.c 
69.3 
609.8 
70.0 
70.0 
70.9 


6.7 
63.9 
75.8 
78.7 
(9,7 
80.0 
80.0 
80.0 
80.0 


awe 
68.6 
83.4 
87.6 
89.0 
ooo 
OF. 9 
90.0 
90.0 
90.0 


8.0 
Plea) 
lee? 
97.8 
99.4 
797 


6.7 
eae 
eh. 9 
147.5 
S30 
209.6 
240.0 
270.0 
S005 9 
530.0 


6.7 
85.0 
We. 
171.9 
207.8 
244.8 
280.0 
315.0 
a0. 0 


Cet 
99.8 
156 6 
197.7 
239.6 
280.0 
320.0 
560.0 
400.0 


Je 
LO Lee 
109.9 
ee0c8 
267.4 
314.5 
300.0 
405.0 
450.0 
495.0 


8.0 
122.5 
186.0 
246.4 
239 
349.0 





400 
450 
500 
550 


ET i i i i a a a a a a a a a, a ee, a ee a a ee. a Se 


OUTPUT FROM DIVIDE AND COLOR PROGRAM 


ee a 
ee 
ooo - 


99.9 
100.0 
100.0 
100.0 


339 wc 
450.0 
500.0 
550.0 


x 


i ee ee ee a a a ee ee, ee, a, ee, ee, ee, ee ee a, ee, a, ee ee, ee, a, a Se | 


N = 25, EDGES = 


ORDER BRANCHES 


avg max 

RANOOM 6 6 
DEGREE 6 6 
DEG 3 6 6 
ADJCEN 6 6 
LCA 6 6 


N = 25, EDGES = 


ORDER BRANCHES 
avg max 

RANDOM 12 l2 
DEGREE 12 12 
DEG 3 12 le 
ADJCEN te 12 
LCA te 12 


a) = o5, EDGES = 
ORDER BRANCHES 


RANDOM 17 18 
MEGREE 17 19 
wes 3 17 18 
ADJCEN 17 18 

LCA 17 18 


N = 25, EDGES = 
ORDER BRANCHES 


avg max 
RANDOM 21 2S 


6( 2%) 
NODES 
avg nax 
12 12 
le le 
le 12 
le 12 
le 
12€ 4%) 
NODES 
avG max 
23 c4 
a) 24 
23 24 
23 e4 
23 24 
18( 6%) 
NODES 
avg Max 
34 56 
34 36 
34 36 
34 36 
34 36 
24( 8%) 
NODES 
avg max 
43 So 





DEGREE ee 46 
DEG 3 ei e7 
ADJCEN 21 23 
BCA e211 23 

N = 25, EDGES = 
ORDER BRANCHES 
avgQ max 

RANDOM 25 4S 
DEGREE e7 52 
DEG 3 24 31 
ADJCEN 23 30 
LCA 23 30 

N = 25, EDGES = 
ORDER BRANCHES 
avg max 

RANDOM Si 945 
DEGREE 34 i9% 
Beco 3 26 94 
ADJCEN 24 35 
LCA 24 35 

N = 25, EDGES = 
ORDER BRANCHES 
avg max 

RANDOM 84 1012 
DEGREE 43 208 
DEG 3 29 86 
ADJCEN 29 170 
meA 25 S4 


46 90 
a4 55 
43 54 
44 54 
30(10%) 
NODES 
avg max 
55 Gi 
Sal Oy 
Se 68 
ae 67 
S2 68 
36(12%) 
NODES 
avg max 
108 1893 
72 380 
S59 189 
57 82 
S57 82 
42(14%) 
NODES 
avQ max 
172 2021 
91 425 
64 177 
65 349 
58 109 


ork xX ke kK kK Ke Ke Ke Ke ke Ke Ke Re Ke eR Ke RR Re Re Re eR he he Re 


RESULTS OF VARYING EDGES FOR A CONSTANT N 


nme ke © eK kK kK kK KR eR RR eh hU Kr hh hl hl hl hlUh hl hUlUr Thr KY RH RR eR 


N= 20, EDGES = 38(20%) 
order branches nodes Logarithm of max 
avQ max avqQ max branches nodes 
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DEGREE 46 363 96 731 S09 
DVEC 3 e7 137 60 278 4.92 
ADJCEN 23 82 Se 166 4.41 
LCA 20 33 46 70 3.50 
mmene0, EDGES = 47(25%) 
order branches nodes 
avg max avqQ max branches 
DEGREE 4&5 t!e7 94 248 4.84 
DVEC 3 28 93 62 18e 4.53 
ADJCEN 21 43 Se 87 554756 
LCA 20 el 48 63 3.30 
N = 20, EDGES = 57(30%) 
order branches nodes 
avgG max avg max branches 
DEGREE 49 187 99 380 Sce5 
DVEC 3 32 164 67 320 Sa) 0 
ADJCEN 25 109 56 223 4.69 
LCA e21 $5 48 117 4.01 
fees €0, EDGES = 66(352%) 
order oranches nodes 
avg max avg max branches 
DEGREE tiie 1119 219 2198 Tae 
OyvetC 5 Se 308 104 608 Si 3 
ADJCEN 38 513 Sompel> 6.e4 
LCA 23 6e 51 123 4.13 
N = 20, EDGES = 76(40%) 
order branches nodes 
avg max avg max branches 
DEGREE NOT ENOUGH SPACE AVAIL 
OVEC 3 59 267 116 524 S99 
ADJCEN 30 292 65 593 5.68 
LCA 23 285 Se 568 a0) 
N = 20, EDGES = 85(45%) 


Se 27 
aor 
5.11 
4.25 


Logarithm of max 


nodes 
Soe 
S.e0 
4.47 
4.14 


Logarithm of max 


nodes 
5.94 
Set 
5.41 
4.76 


Logarithm of max 


nodes 
for C 
6.41 
6.92 
4.81 


Logarithm of max 


nodes 


6.26 
6.39 
6.34 





order 


DEGREE 
OVEC 3 
ADJCEN 

LCA 


N = 20, 


order 


DEGREE 
DVEC 3 
ADJCEN 

ECA 


Ns 
order 


DEGREE 
DVEC 3 
ADJCEN 

LCA 


N = 20, 


order 


DEGREE 
BVEC $ 
ADJCEN 

LCA 


N = 20, 


order 


DEGREE 
DVEC 3 
ADJCEN 

LCA 


a = 


20, 


20, 


oranches 
avg max 
Ge 635 
So 425 
37 06379 
23 64 
EDGES = 
branches 
avg max 
101 B11 
S7 367 
34 S01 
e8 584 
EDGES = 
oranches 
avg max 
106 431 
69 440 
34 258 
23 79 
EDGES = 
branches 
avgQ max 
99 969 
62 576 
37 =6561 
ef 438 
EDGES = 
branches 
avg max 
86 765 
S55 468 
28 188 
el 53 
EDGES = 


nodes 
avqQ max 
173 1206 
109 818 
77 06749 
49 128 


95(50%) 


nodes 
avg max 
188 1608 
107 676 
71 £4980 
60 1150 


104(55%) 


nodes 
avg max 
194 777 
131 844 
70 497 
48 156 


114(60%) 


nodes 
aVG max 
180 1696 
114 1029 
75 1088 
S56 855 


123(65%) 


nodes 
avG max 
155 1351 
101 83e 
S57 365 
43 103 


PSS (70%) 


Logarithm of max 


oranches nodes 
6.45 7.10 
6.05 Ga tf 
5.94 6.62 
4.16 4.85 


Logarithm of max 


branches nodes 
6.70 aa 
S91 6.52 
Seee Go 
60 Lets 


Logarithm of max 


branches nodes 
6.07 6.66 
6.09 6.74 
555 6.21 
4.37 S205 


Logarithm of max 


branches nodes 
S.1o0 7.44 
6.36 6.94 
6.33 0.99 
6.08 Oat 


Logarithm of max 


branches nodes 
6.64 pee 
6.15 6.72 
5.24 Se 0 
3.97 4.63 
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order branches 
avg max 

DEGREE 74 419 
DVEC 3 47 282 
ADJCEN 21 57 
LCA 20 55 

N = 20, EDGES = 
order branches 
avg max 

DEGREE 60 295 
OVEC 3 38 226 
ADJCEN 2e 66 
LCA 19 37 

N = 20, EDGES = 
order branches 
avg max 

MeGREE 45 143 
DVEC 3 30 at 
ADJCEN 20 46 
LCA 19 eT 


nodes 
avg max 
130 748 
86 498 
GG 114 
4} 62 
142(75%) 
nodes 
avqQ max 
103 499 
69 403 
43 118 
39 7i 
152(80%) 
nodes 
avg max 
76 237 
sa) 118 
39 81 
S71 4g 


Logarithm of max 


branches 
6.04 
5.64 
4.04 
3.50 


nodes 
6.6¢ 
6.21 
Q.74 
4.13 


Logarithm of max 


branches 
$5.69 
5.42 
4.19 
3.61 


nodes 
Oe | 
6.00 
a7 7 
4.26 


Logarithm of max 


branches 
4.96 
4.26 
Saa5 
5.230 


Wet 


nodes 
Se 
Qi/77 
4.39 
559 
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